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CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

The  pavement  structure  is  the  single  most  costly  element  of  America’s  highway 
system  (AASHTO  Guide,  1993).  It  is  not  surprising  then  that  Federal,  State,  and  local 
Transportation  Departments  are  continuously  seeking  advances  in  technology  to  reduce 
this  cost.  The  development  of  mechanistic-empirical  approaches  for  rehabilitation 
designs  based  on  layered  elastic  theory  is  very  promising.  The  falling  weight 
deflectometer  (FWD)  is  the  device  most  often  used  to  provide  the  field  data  necessary  for 
mechanistic-empirical  calculations.  As  mechanistic-empirical  approaches  and  layered 
elastic  theory  design  aids  become  more  commonly  used,  pavement  design  will  steadily 
transition  from  an  art  to  a  science  based  on  accepted  theory. 

The  ultimate  goal  is  to  accurately  depict  the  remaining  life  of  a  given  pavement 
section  and  to  design  an  overlay  that  is  neither  over-conservative  nor  under-conservative. 
These  analysis  methods  can  also  be  used  to  prioritize  pavements  for  rehabilitation  based 
on  scientific  data  as  opposed  to  general  rules-of-thumb.  The  pavement  engineer  will  then 
have  a  clear  picture  of  the  priority  for  pavement  overlays  and  an  accurate  measure  of  the 
thickness  required  to  extend  the  life  of  the  pavement. 

Mechanistic  design  procedures  are  based  on  layered-elastic  theory  which 
supposes  that  pavements  can  be  modeled  as  a  multi-layered  elastic  structure  placed  on  an 
elastic  foundation.  If  pavements  are  modeled  in  this  manner  it  is  possible  to  calculate 
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stresses,  strains  and  deflections  due  to  traffic  and/or  environmental  conditions  at  any 
point  within  or  below  the  structure  (AASHTO  Guide,  1993).  Because  pavement 
performance  is  also  affected  by  other  parameters  which  cannot  be  directly  modeled  by 
mechanistic  models,  empirical  correlations  must  be  incorporated  into  the  process.  This  is 
why  it  is  referred  to  as  a  mechanistic-empirical  design  procedure. 

1.2  SCOPE 

The  goal  of  the  research  described  in  this  thesis  is  to  develop  relationships  using 
measurable  pavement  responses  to  loading  to  predict  the  current  pavement  life  and  the 
overlay  requirements.  The  equations  are  designed  to  be  used  with  any  pavement 
nondestructive  testing  device  that  can  measure  the  pavement  deflection  basin.  The  most 
popular  device  currently  in  use  is  the  FWD.  The  most  promising  new  technology  is  a 
laser  rolling  wheel  deflectometer;  prototypes  are  being  developed  by  Quest  Integrated 
and  the  Swedish  National  Road  Administration. 

1.3  STRUCTURE 

This  thesis  is  structured  as  follows: 

Chapter  1  is  the  Introduction. 

Chapter  2  is  the  Literature  Review.  Its  purpose  is  to  describe  previous  research 
conducted,  basic  theories  and  applications,  and  current  practices  pertinent  to  this  thesis. 

Chapter  3  describes  the  research  conducted  and  the  methodology  for  developing 
predictive  relationships  for  life  remaining  in  the  pavement  structure. 
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Chapter  4  describes  the  research  conducted  and  the  methodology  for  developing 
an  adjustment  factor  for  load  to  adjust  both  Area  Parameter  and  Dq  to  a  standard  load. 

Chapter  5  describes  the  research  conducted  and  the  methodology  for  developing 
predictive  relationships  for  the  overlay  required  on  a  pavement  system. 

Chapter  6  provides  some  conclusions  and  recommendations  for  further  research 
in  this  area. 

The  Appendices  contain  much  of  the  data,  charts,  and  equations  that  are 
important  background  to  this  thesis  but  were  not  necessary  in  the  main  text. 


CHAPTER  2 
LITERATURE  REVIEW 

2.1  INTRODUCTION 

This  chapter  consolidates  background  information  on  nondestructive  testing  and 
design  procedures  using  nondestructive  testing  results.  Research  efforts  similar  to  that 
described  in  this  thesis  are  also  presented. 

2.2  NONDESTRUCTIVE  TESTING 

Nondestructive  testing  (NDT)  is  a  term  used  to  describe  many  different 
techniques  to  analyze  the  properties  of  a  material  or  structure  without  altering  the  in-situ 
properties  of  the  material  or  structure  in  any  way.  For  obvious  reasons  NDT  is  the 
preferred  method  to  test  a  wide  variety  of  structures  from  high-pressure  steam  lines  in 
nuclear  powered  submarines  to  pavements.  NDT  is  a  growth  area  for  many  emerging 
technologies  as  greater  accuracy  and  precision  are  required  in  a  wide  variety  of 
industries.  For  the  remainder  of  this  thesis,  NDT  will  refer  only  to  the  analysis  of 
pavement  structures,  and  more  specifically  to  pavement  surface  deflection  measurement 
devices. 

The  most  reliable  method  for  determining  the  material  characteristics  in  an 
existing  pavement  system  is  to  actually  retrieve  a  sample  by  coring  through  every  layer 
and  conducting  laboratory  tests  on  each  layer,  but  costs  and  time  are  a  problem.  In 
addition,  samples  of  unstabilized  material  are  very  difficult  to  collect  and  preserve  so  in- 
situ  properties  are  unaffected.  For  these  reasons,  different  methods  of  nondestructive 
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testing  have  been  developed  in  an  effort  to  maximize  the  accuracy  and  reliability  of  the 
data  while  minimizing  the  cost  to  conduct  the  test. 

Current  NDT  devices  for  analyzing  pavement  structures  can  be  divided  according 
to  the  types  of  loading  placed  on  the  pavement.  Static  deflection  devices  like  the 
Benkleman  Beam,  the  Plate  Bearing  Test  and  the  LaCroix  Deflectograph  were  the 
pioneering  efforts  in  this  field.  Steady  state  deflection  equipment  (Dynaflect  and  Road 
Rater  Models  400  and  2000)  uses  a  vibratory  load  and  is  still  used  by  many  agencies 
(Mahoney  et  al.  1995;  WSDOT  Pavement  Guide,  1995).  Current  state  of  the  art 
technology  utilizes  an  impact  (or  impulse)  load  which  approximates  a  vehicular  load  on 
the  pavement  better  than  a  static  or  vibratory  load  (ASTM  D4694-87,  1989).  These 
impact  load  devices  are  commonly  referred  to  as  Falling  Weight  Deflectometers  (FWDs) 
and  are  manufactured  by  several  companies  worldwide  including  KUAB,  Dynatest, 
Foundation  Mechanics,  and  Phoenix.  Over  the  last  decade,  the  FWD  has  become  the 
deflectometer  of  choice  for  most  road  authorities  worldwide;  hence  there  has  been 
considerable  analysis  of  FWD  data  (Mahoney,  et  al.,  1995).  Table  2-1  shows  the  past, 
current  and  predicted  usage  of  the  various  NDT  devices  currently  used  in  the  United 
States.  Appendix  E  provides  an  overview  of  several  of  the  most  common  nondestructive 
testing  devices;  past,  present,  and  future.  Further  information  is  available  in  a  synthesis 
study  of  NDT  devices  by  Smith  and  Lytton  (1984).. 
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Table  2-1  State  Highway  Administration  (SHA)  Deflection  Equipment 
(WSDOT  Pavement  Guide,  1995) 


Device 

Number  of  SHAs  Using  Device  for 
Various  Time  Periods 

mid-1980s 

1990 

mid-1990s 

(estimated) 

Benkelman  Beam 

18 

3 

0-3 

Dynaflect 

18 

11 

5-10 

Road  Rater 

5 

4 

2-6 

FWD 

5 

30 

30-40 

2.3  INTERPRETATION  OF  NDT  DATA 

NDT  data  is  typically  converted  into  different  deflection  basin  parameters  for  use 
in  analyzing  the  pavement.  A  summary  of  deflection  basin  parameters  is  found  in  the 
WSDOT  Pavement  Guide  for  Design  Evaluation,  and  Rehabilitation  (Feb.  1995)  and  is 
shown  in  Table  2-2. 

Two  important  deflection  basin  parameters  are  Dq  and  Area  Parameter.  Because 
these  are  used  extensively  throughout  this  thesis,  special  attention  will  be  paid  to  them 
here.  Dq  is  the  maximum  deflection  measured  at  the  center  of  the  load  applied  to  the 
pavement.  It  is  measured  in  units  of  length.  The  Area  Parameter,  A,  is  a  measurement 
of  the  normalized  area  of  a  slice  of  the  deflection  basin  from  the  center  of  the  load  to 
90.5  cm  (3  ft)  from  the  center  of  the  load.  Because  it  is  normalized  by  the  Dq  value,  it 
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has  units  of  length.  It  is  calculated  using  the  equation  in  Table  2-2,  using  the  maximum 
deflections  for  Dq,  Dj,  Dj,  and  D3.  Figure  2-1  gives  a  graphical  presentation  of  how 
Area  Parameter  is  measured  and  what  it  actually  represents.  The  maximum  value  to 
expect  for  Area  Parameter  is  915  mm  (36.0  in)  which  would  indicate  all  deflections  were 
equal.  This  would  represent  an  extremely  stiff  pavement  system.  The  minimum  value  to 
expect  for  Area  Parameter  is  280  mm  (11,1  in).  The  lower  range  of  Area  Parameter 
values  indicate  that  the  pavement  is  not  acting  significantly  differently  than  the 
underlying  layers  (WSDOT  Pavement  Guide,  1995).  Table  2-3  provides  a  guide  to 
pavement  type  and  condition  typically  associated  with  Area  Parameter  values. 

Deflection  basin  parameters  can  be  used  in  several  sets  of  equations  (shown  in 
Table  2-3)  to  compute  the  subgrade  elastic  modulus.  The  accuracy  of  these  methods  is 
far  better  for  the  subgrade  than  for  the  surfacing  layer  (WSDOT  Pavement  Guide,  1995). 
Surfacing  layer  equations  can  be  found  in  the  WSDOT  Pavement  Guide  (1995). 


Table  2-2  Summary  of  Deflection  Basin  Parameters  (WSDOT  Gnide,  1995) 


Parameter 

Formula 

Measuring  Device 

Maximum  Deflection 

Do 

Benkelman  Beam 

LaCroix  Deflectograph 

FWD 

Radius  of  Curvature 

R  =  rV{2Do(Do/D,-l)} 

r=5 

Curvaturemeter 

Spreadability 

S  =  100[(Do+Di=D2+D3)/5]/Do 

D1...D3  spaced  12  apart 

Dynaflect 

Area 

A  =  6[1  +  2(Di/Do)  +  2(D2/Do) 

+  D3/D0)];  0,1,2,3  ft 

FWD 

Shape  Factors 

Fi  ~  (Do-D2)/Di 

F2  =  (Dj-D3)/D2 

FWD 

Surface  Curvature  Index 

SCI  =  Do-Dp  where 

r  =  12  or  r  =  20 

Benkelman  Beam 

Road  Rater 

FWD 

Base  Curvature  Index 

BCI  =  D24»"  D36" 

Road  Rater 

Base  Damage  Index 

BDI  ~  “  D24" 

Road  Rater 

Deflection  Ratio 

Qj.  =  D/Do,  where  Dr~Do/2 

FWD 

Bending  Index 

BI  =  D/a,  where 

a  =  Deflection  basin 

Benkelman  Beam 

Slope  of  Deflection 

SD  =  tan'^(Do-Dr)/r 
r  =  24" 

Benkelman  Beam 
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0  305mnn  610mm  915mm 


A 


Equal  Area  Bounded  by 
Do  and  Area  Parameter 


Figure  2-1  Graphical  Representation  of  Area  Parameter 
(after  WSDOT  Pavement  Guide) 


Table  2-3  Area  Parameter  Values  for  Typical  Pavements 
(after  WSDOT  Pavement  Guide) 


Pavement 

Area  Parameter 

mm 

in 

Thick  ACP  (>  107  mm  ACP) 

530-760 

21-30 

Thin  ACP  (<  107  mm  ACP) 

410-530 

16-21 

BST  flexible  pavement  (relatively  thin) 

380-430 

15-17 

Weak  BST 

300-380 

12-15 

Table  2-4  Predictive  Relationships  for  Subgrade  Modulus 


WSDOT  Two  Layer 

Esg=  -466  +  0.00762  (P/D3) 

Esg  =  -198  +  0.00577  (P/D4) 

Esg  =  -371  +  0.00671  (2P/(D3+D4)) 

WSDOT  Three  Layer 

Esg  =  -530  +  0.00877  (P/D3) 

Esg  =  -1 1 1  +  0.00577  (P/D4) 

Esg  =  -346  +  0.00676  (2P/(D3+D4)) 

AASHTO 

Eso'CPXSrVCD.Xr) 

M„  =  P(lV)/(7t)(D,)(r) 

South  African 

logioEsG  “  9727  -  0.989  logio62ooo 

CBR  Correlation 

Esg  (psi)  =  1500(CBR) 

♦Valid  for  CBR  values  less  than  10 

only* 

Esg  (MPa)  =  lO(CBR) 

Definitions  for  variables  used  in  Table  2-4: 


P 

Esg 

Do 

Do.67 

Di 

D2 


=  applied  load  (lbs)  on  1 1 .8  in.  plate 
=  subgrade  modulus  (psi)  (Pa  for  South  African  eqn.) 
=  deflection  beneath  center  of  load  plate  (in.) 

=  deflection  0.67  feet  from  center  of  load  plate  (in.) 

=  deflection  1  foot  from  center  of  load  plate  (in.) 

=  deflection  2  feet  from  center  of  load  plate  (in.) 
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Ds 

D4 

Sf 

D, 

r 

Mr 

^2000 

CBR 


=  deflection  3  feet  from  center  of  load  plate  (in.) 

=  deflection  4  feet  from  center  of  load  plate  (in.) 

=  subgrade  modulus  reduction  factor  based  on  Poisson’s  ratio,  p, 
=  deflection  measured  r  feet  from  the  center  of  load  plate  (in.) 

=  distance  from  the  applied  load  to  (in.) 

=  backcalculated  subgrade  resilient  modulus  (psi) 

=  deflection  2000  mm  from  the  center  of  the  load  plate  (pm) 

=  California  Bearing  Ratio 
=  Poisson’s  Ratio 


2.4  DESIGN  PROCEDURES 


2.4.1  Empirical 

An  empirical  design  approach  concerns  itself  only  with  the  results  observed  in  the 
laboratory  or  in  the  field.  Correlations  are  made  from  large  numbers  of  observations  to 
create  empirical  relationships  between  inputs  and  results.  Hveem  and  Carmany’s 
equation  for  the  thickness  of  asphalt  pavement  cover  required  is  an  example  of  an 
empirical  pavement  design  (WSDOT  Pavement  Guide,  1993).  The  AASHTO  process  is 
another  example  of  empirical  design. 

2.4.2  Mechanistic 

Mechanistic  approaches  to  design  apply  laws  of  physics  and  material  properties  to 
determine  a  structure’s  reaction  to  loading.  Mechanistic  design  procedures  for 
pavements  rely  on  the  assumption  that  a  pavement  can  be  modeled  as  a  multi-layered 
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elastic  or  visco-elastic  structure  on  a  similar  elastic/visco-elastic  foundation.  This 
assumption  allows  the  calculation  of  stresses,  strains,  and  deflections  throughout  the 
pavement  structure  and  subgrade  (AASHTO  Guide,  1993). 

2.4.3  Mechanistic-Empirical 

As  the  name  implies,  a  mechanistic-empirical  approach  incorporates  elements  of 
both  empirical  and  mechanistic  methods.  The  mechanistic  component  uses  mathematical 
models  to  compute  pavement  reactions  of  interest  while  the  empirical  component  relates 
these  reactions  to  some  measure  of  the  performance  of  the  pavement  such  as  remaining 
life  of  the  pavement  (WSDOT  Guide,  1995). 

There  are  many  benefits  to  using  a  mechanistic-empirical  approach  in  pavement 
design.  Because  these  methods  are  based  on  long  established  linear-elastic  theory,  they 
will  model  the  pavement  more  correctly  than  the  empirical  equations  which  have  been 
traditionally  used  for  flexible  pavement  design  (AASHTO  Guide,  1993).  Other  benefits 
(WSDOT  Guide,  1995)  include: 

•  ability  to  accommodate  changing  load  types 

•  better  use  of  available  materials 

•  ability  to  accommodate  new  materials 

•  improved  reliability  of  performance  predictions 

•  better  defined  role  of  construction 

•  improved  definition  of  existing  pavement  layer  properties 

•  ability  to  accommodate  environmental  and  aging  affects  on  materials 


13 


•  material  properties  which  relate  better  to  actual  pavement  behavior  and  performance 

2.4.4  Backcalculation 

Backcalculation  is  a  mechanistic-empirical  evaluation  which  iteratively  computes 
surface  deflection  basins  for  different  moduli  until  the  calculated  basin  matches  a 
measured  deflection  basin  within  a  preprogrammed  tolerance.  The  layer  moduli  that 
produce  the  calculated  deflection  basin  can  be  limited  to  certain  ranges  to  ensure  they 
make  engineering  sense  (WSDOT  Pavement  Guide,  1995).  Backcalculation  is  achieved 
by  using  different  software  programs  designed  to  be  used  on  standard  microcomputers 
due  to  its  calculation  intensity.  These  programs  usually  use  a  linear-elastic  analysis. 

A  successful  analysis  of  a  pavement  based  on  a  measured  deflection  basin 
depends  on  the  accuracy  of  several  input  parameters,  including  the  layer  thicknesses  and 
composition  (which  will  dictate  the  initial  moduli  and  Poisson’s  Ratio  to  be  used) 
(Lytton,  et  al.,  1990).  When  these  parameters  cannot  be  determined  from  construction 
records,  core  samples  should  be  drilled.  Some  researchers  and  agencies  have  had  success 
using  seismic  refraction  studies  or  ground  penetrating  radar  for  determining  layer 
thicknesses. 

2.4.5  Computer  Pro2rams 

The  Federal  Aviation  Administration  (FAA)  has  developed  a  program  called 
LEDF AA  which  makes  some  standard  assumptions  concerning  the  materials  and  the 
structural  model  of  the  pavement  system  (FAA  Workshop,  1995).  These  assumptions. 
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detailed  below,  are  consistent  with  most  other  software  programs  using  layered  elastic 
theory: 

•  Multi-layer  elastic  system 

•  All  layers  are  infinite  in  the  lateral  direction 

•  Each  layer  except  the  bottom  layer  has  a  constant,  finite  thickness 

•  Friction  at  the  surface  is  not  considered  (surface  shearing  neglected) 

•  The  bottom  layer  is  infinite  in  depth 

•  Each  layer  consists  of  a  homogenous  and  isotropic  material 

•  The  elastic  modulus,  E,  and  Poisson’s  ratio,  V,  for  each  layer  are  used  to  calculate  the 

deflection  basin  (pavement  response  to  loading) 

The  Washington  State  Department  of  Transportation  has  developed  three 
computer  programs  to  enhance  their  pavement  management  and  design  capability. 
EVERSTRS,  a  layered  elastic  analysis  program,  is  used  to  estimate  the  stresses,  strains 
and  deflections  within  a  layered  pavement  system  due  to  a  static  load.  EVERCALC,  a 
pavement  analysis  computer  program,  is  used  to  estimate  the  “elastic”  moduli  of 
pavement  layers.  Given  the  FWD  load,  sensor  spacing,  seed  moduli  (initial  guesses)  and 
layer  thicknesses,  EVERCALC  computes  an  estimate  of  the  actual  moduli  by  iteratively 
comparing  the  computed  deflection  basin  with  the  measured  deflection  basin  until  they 
match  within  given  tolerance  parameters.  Figure  2-2  shows  a  simplified  flow  chart  for 
EVERCALC.  EVERPAVE,  an  empirical-mechanistic  overlay  design  program,  is  used  to 
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estimate  the  overlay  thickness  required  for  a  pavement  based  on  fatigue  cracking  criteria 
and  rutting  failure  criteria.  These  programs  are  extremely  useful  tools,  but  should  not 
replace  the  experience  and  judgment  of  a  pavement  engineer  to  verify  the  reasonability 
of  the  results. 


Figure  2-2  Simplified  Flow  Chart  for  EVERCALC 
(after  WSDOT  Pavement  Guide, 1995) 


2.4.6  Regression  Equations 

Regression  equations  (or  algorithms)  can  be  developed  for  data  that  correlates 
well  enough  to  provide  “predictive”  relationships.  Other  researchers  have  studied  and 
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created  regression  equations  for  the  same  or  similar  parameters  that  this  thesis  deals  with. 
Some  background  information  about  regression  and  the  meaning  of  some  of  the 

descriptors  such  as  and  Root  Mean  Square  Error  (RMSE)  are  discussed  in  Section 
3.2.2  and  Appendix  C. 

Gomez-Achecar  and  Thompson  (1986)  developed  regression  equations  relating 
pavement  characteristics  and  deflections  to  certain  strains,  stresses  and  deflections.  They 
are  presented  in  Table  2-5.  In  general,  the  AC  thickness  and  modulus  must  be  known  in 
addition  to  the  subgrade  modulus  to  use  these  equations.  The  R^  and  RMSE  are  very 
good  for  the  equations  using  the  three  independent  variables  just  mentioned.  The  goal  of 
this  thesis  IS  to  create  very  similar  relationships  that  do  not  require  knowledge  of  the  AC 
modulus  or  thickness. 

Olle  Andersson,  Professor  of  Highway  Engineering  at  the  Royal  Institute  of 
Technology  in  Sweden,  presented  data  from  140  different  pavement  scenarios  (varying 
AC  thickness,  AC  Modulus,  Base  Thickness,  Base  Modulus,  Subgrade  Modulus).  His 
findings  (Andersson,  undated  paper)  include  graphs  indicating  a  linear  relationship 
between  the  displacement  at  the  center  of  the  load  and  other  displacement  parameters  for 

given  AC  modulus  s  and  subgrade  modulus’s.  No  regression  equations  were  proposed  in 
this  paper. 
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Table  2-5  Regression  Equations  Developed 
by  Gomez-Achecar  and  Thompson  (1986) 


Equation 

■ 

RMSE 

Asphalt  Strain 

Log  Sac  ~  5.746  -  1.589  Log  T^c  -  0.774  Log  E^c  -  0.097  Log  Ejy 

0.967 

0.083 

Suberade  Deviator  Stress 

Log  DEV  =  2.744  -  1,138  Log  Tac  -  0.515  Log  Eac  +  0.289  Log  Er; 

0.976 

0.053 

Surface  Deflection 

Log  Do  =  3.135  -  0.895  Log  Tac  -  0.359  Log  Eac  -  0.287  Log  Erj 

0.984 

0.033 

Suberade  Deflection 

Log  DSUB  =  3.090  -  0.979  Log  Tac  -  0.321  Log  Eac  -  0.306  Log 

Eri 

0.983 

0.037 

Suberade  Vertical  Strain 

Log  s,  =  4.022  -  1.680  Log  Tac  -  0.667  Log  Eac  -  0. 165  Log  Eri 

0.944 

0.110 

Asphalt  Strain- Siuf ace  Deflection 

Log  Sac  =  1.53  Log  Do  +  0.3 19 

0.81 

0.20 

Suberade  Stress  Ratio-Surface  Deflection 

none 

none 

Log  S=  1.28  Log  Do -2.21 

given 

given 

where: 

Tac  ~ 

Asphalt  Concrete  Thickness  (inches) 

Eac  “ 

Asphalt  Concrete  Modulus  (ksi) 

Eri  = 

Subgrade  Modulus  (ksi) 

Sac  “ 

Asphalt  Concrete  Radial  Tensile  Strain  (microstrain) 

DEV  = 

Subgrade  Deviator  Stress,  Oj  -  03  (psi) 

Do  = 

Surface  Deflection  (mils) 
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DSUB  =  Subgrade  Deflection  (mils) 

Sj,  =  Subgrade  Vertical  Strain  ( X 

S  =  Subgrade  Stress  Ratio  (subgrade  deviator  stress/unconfmed  compressive 

strength) 

2.5  SOURCES  OF  ERROR 

There  are  two  primary  sources  of  error  in  the  use  of  NDT  data  to  compute  layer 
moduli.  These  are  random  errors  and  systematic  errors  which  are  discussed  in  this 
section. 

2.5.1  Random  Errors 

Random  errors  are  typically  errors  in  the  measurements  themselves  either  due  to 
the  measurement  sensors  or  the  load  cells.  These  errors  cannot  easily  be  predicted  but 
they  can  be  reduced  by  repeating  measurements  and  averaging  the  results  (Lytton,  et  al. 
1990).  For  Laser  rolling  wheel  deflectometers,  random  errors  should  not  affect  the 
results  because  of  the  vast  amount  of  measurements  which  will  be  averaged  over  a  small 
section  of  pavement. 

2.5.2  Systematic  Errors 

Systematic  errors  are  by  far  the  most  troublesome  because  they  are  more  difficult 
to  correct,  are  numerous,  and  have  the  tendency  to  compound  their  effects.  Systematic 
errors  are  the  result  of  erroneous  assumptions  made  in  the  backcalculation  process. 
Lytton,  et  al  (1990)  described  research  which  included  the  creation  of  an  ‘expert  system’ 
to  be  used  in  conjunction  with  a  backcalculation  computer  program. 
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An  engineer’s  knowledge  of  pavements  is  critical  to  the  mechanistic-empirical 
process.  This  process  cannot  be  used  as  a  ‘cookbook’  to  generate  solutions  without 
proper  analysis  of  the  results.  The  AASHTO  Guide  (1993)  notes  that  good  engineering 
judgment  on  the  part  of  the  designer  will  always  be  critical,  especially  since  so  many  of 
the  inputs  to  the  design  cannot  be  modeled  properly  by  mechanistic  models. 

2.6  SIGNIFICANT  PARAMETERS  AFFECTING  NDT 

Pavement  systems  react  very  differently  to  different  loading  conditions.  The 
magnitude  and  the  rate  of  loading  can  dramatically  alter  the  pavement  response 
(deflection).  Since  most  NDT  measures  the  pavement’s  response  to  loading  in  order  to 
estimate  the  pavement  strength  properties,  it  is  extremely  important  to  keep  the  loading 
conditions  consistent. 

2.6.1  Stress  Sensitivity  of  Pavement  Materials 

The  most  significant  factors  that  affect  the  response  of  the  asphalt-concrete 
layer(s)  are  the  rate  or  frequency  of  loading  and  the  temperature  (Hoffman  and 
Thompson,  1981).  The  AC  layer  is  not  greatly  influenced  by  the  stress  imposed  on  it. 

The  unbound  materials  often  found  in  the  base  course,  subbase  and  subgrade  have 
been  found  to  be  sensitive  to  the  stress  exerted  upon  them.  Generally,  the  modulus  of 
fine-grained  materials  decrease  with  increasing  stresses  while  the  modulus  of  granular 
materials  increase  with  increasing  applied  stress  (Newcomb,  1986).  Since  the 
unstabilized  base  course,  subbase  and/or  subgrade  can  respond  in  a  nonlinear  fashion 
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with  regards  to  applied  stress,  the  modulus  may  be  calculated  by  either  of  two  equations 
(WSDOT  Guide,  1995): 

E  =  K,(0)^ 
or 

E  =  K3(aa)^ 

where: 

E  =  Layer  Modulus 

Ki  K4  =  Regression  Constants 

6  =  Sum  of  the  principal  stresses  (oj  +  02  +  O3) 

Od  =  Difference  between  major  and  minor  principal  stresses  (Oi  -  03) 

The  laboratory  method  typically  used  to  compute  the  regression  constants  is  a  repeated 
load  triaxial  test  on  disturbed  samples  which  are  recompacted  (Newcomb,  1986; 
Hoffman  and  Thompson,  1981). 

Studies  conducted  on  lime  and  cement  treated  materials  (Hoffman  and 
Thompson,  1981)  indicate  that  their  resilient  response  lies  somewhere  between  a 
completely  bonded  material  (elastic)  and  the  unbound  raw  material  (nonlinear)  and 
depend  a  great  deal  on  factors  including  mix  proportions,  curing  time,  method  of  testing 


and  others. 
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2.6.2  Variations  in  Loading  Type  and  Duration 

Three  different  mathematical  formulations  can  be  used  to  determine  the  response 
of  a  pavement  system  under  different  loading  conditions.  The  are  (Hoffman  and 
Thompson,  1986); 

1 .  Layered  Theory:  Does  not  consider  time  or  mass 

2.  Visco-elastic:  Considers  time  only 

3.  Dynamic:  Considers  time  and  mass 

One  model  should  not  be  used  to  predict  results  for  another.  They  are  different  models 
and  are  not  interchangeable.  The  research  described  in  this  thesis  used  the  layered 
theory. 

During  the  1950s  and  1960s  researchers  noted  that  pavement  surface  deflections 
decrease  as  vehicle  speeds  increase  (Hoffman  and  Thompson,  1981).  Laser  Rolling 
Wheel  Deflectometers  will  experience  this  effect;  consequently,  dynamic  analysis 
methods  should  be  developed  to  properly  interpret  the  data  generated  by  these  machines. 

2.6.3  Temperature  Effects 

Because  asphalt  properties  are  highly  temperture-dependent,  the  temperature  of 
the  asphalt  layer  affects  the  deflections  measured  by  any  NDT  device.  Many  researchers 
have  attempted  to  explain  and  predict  the  relationship  between  temperature  and 
deflection.  It  is  very  important  to  adjust  deflections  to  a  standard  temperature  before 
attempting  to  use  the  deflections  in  any  kind  of  established  analysis  or  design.  Previous 
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research  has  yielded  several  important  general  conclusions,  some  of  which  are  outlined 
here. 

The  Canadian  Good  Roads  Association  (Sebastyan,  1961)  determined  that: 

1.  Temperature  variations  could  not  be  correlated  with  mix  design  characteristics.  A 
linear  relationship  was  assumed  between  temperature  and  deflection.  For  this  linear 
relationship,  a  slope  of  0.002  inches  per  10  degrees  Fahrenheit  (pavement  surface 
temperature)  worked  well  for  a  variety  of  pavements. 

2.  No  temperature  effects  could  be  seen  on  pavement  systems  with  very  thick  bases. 

Welch  (1961)  further  concluded  that  the  “deflection  correction  factor  in  inches 
for  a  weak  base  was  considerably  higher  than  that  for  a  stronger  base”. 

During  the  AASHO  Road  Test  (Benkelman,  et  al,  1962)  the  following 
conclusions  were  made; 

1 .  Temperature  has  the  greatest  effect  on  deflection  in  lower  temperature  ranges. 

2.  Very  little  temperature  effect  is  experienced  when  the  average  pavement  temperature 
is  above  26.7  ^  C  (80  °  F). 

In  addition  to  conclusions  made  based  on  observations  of  test  sites,  Meyerhoff  proposed 
the  following  mathematical  relationship  between  deflection  and  temperature  based  on  a 
two-layer  elastic  system  (Meyerhoff,  1962): 

.  0.52W 


where  d 


deflection  in  inches 
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W  =  Wheel  load  in  pounds 

Eg  =  Subgrade  Elastic  Modulus 

Ep  =  Pavement  Elastic  Modulus 

t  =  Pavement  Thickness  in  inches 

Kingham  (1969),  summarized  the  significant  findings  from  many  research  efforts  as 
follows: 

1 .  Surface  pavement  temperatures  do  not  correlate  well  to  deflection  for  asphalt  layers 
greater  than  3  inches.  Mean  temperature  throughout  the  slab  provides  a  better  correlation 
for  these  thicker  pavements. 

2.  Pavements  with  high  deflections  experience  greater  temperature  effects  (thin 
pavements,  weak  subgrade,  etc.). 

3.  Temperature  vs.  deflection  relationships  are  curvilinear  for  large  temperature  ranges 
(greater  than  4.4  °  C  (40  °  F)).  For  smaller  temperature  ranges,  such  as  those  experienced 
in  a  24  hour  period,  the  relationship  can  be  reasonably  modeled  as  a  linear  relationship. 

4.  Field  measurements  do  not  fully  support  theoretical  models  for  temperature  effects. 
Theory  should  not  be  the  only  input  to  a  temperature  correction  procedure. 

Kingham  went  on  to  outline  a  temperature  conversion  procedure  which  can  be  used  to 
normalize  measured  deflections  to  21 . 1  °  C  (70  °  F). 

The  Asphalt  Institutes’  Manual  Series  No.  17,  “Asphalt  Overlays  for  Highway 
and  Street  Rehabilitation”,  (1983)  included  a  discussion  of  adjusting  a  Representative 
Rebound  Deflection  (RRD)  to  a  standard  temperature.  The  RRD  was  based  on 
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Benkelman  Beam  results.  The  manual  provides  a  figure  which  provides  a  temperature 
adjustment  factor  for  Full-Depth  and  three-layered  asphalt  concrete  pavements. 

Thickness  of  the  untreated  aggregate  base  and  mean  pavement  temperature  must  be 
known  to  use  this  chart. 

A  temperature  adjustment  can  be  made  to  Asphalt  Concrete  Modulus,  E^c  with 
the  following  equation  (WSDOT  Pavement  Guide,  1995); 

„  _  9.1242  -  0.04456(T) 

where  E^c  =  AC  Modulus  (MPa) 

T  =  mean  asphalt  concrete  temperature  (degrees  Celsius) 

2.7  DEVELOPING  STANDARD  PROCEDURES 

In  1986,  Harold  Von  Quintas  presented  a  first  draft  of  a  “Standard  Practice  for 
Calculating  In-Situ  Resilient  Modulus  of  Pavement  Materials”  in  an  attempt  to  develop  a 
standard  for  a  procedure  that  is  gaining  acceptance  and  continues  to  evolve  (May  and 
Von  Quintas,  1994).  Seven  drafts  and  six  years  later  the  work  on  the  proposed  standard 
was  put  on  hold  until  after  the  1993  Symposium  “Nondestructive  Testing  of  Pavements 
and  Backcalculation  of  Moduli”  sponsored  by  the  American  Society  for  Testing  and 
Materials  (ASTM).  To  date,  there  seems  to  be  very  little  agreement  on  when  and  how  to 
standardize  backcalculation  procedures. 
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2.8  METHODS  EMPLOYED  BY  STATE  HIGHWAY  ADMINISTRATIONS  AND 
OTHER  AGENCIES 

In  1991  a  National  Cooperative  Highway  Research  Program  (NCHRP)  survey 
was  conducted  to  quantify  state  and  Canadian  Province  standard  practices  in  design 
(WSDOT  Guide,  1995).  Figures  2-3  to  2-6  show  pertinent  results  concerning  methods 
used  by  SHAs  and  other  agencies  to  determine  different  pavement  parameters.  The  data 
to  create  these  charts  came  from  the  WSDOT  Pavement  Guide,  Volume  2,  1995. 


Figure  2-3  Agencies  Use  of  Design  Procedures 


Percentages  are  based  on  63  agencies  reporting. 
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Figure  2-4  Agencies  Method  of  Asphalt  Concrete  Characterization 

Percentages  based  on  51  agencies,  updated  to  include  WSDOT. 


Elastic  Modulus 
11% 


52% 


Figure  2-5  Agencies  Method  of  Base/Subbase  Characterization 

Percentages  based  on  27  agencies  reporting  use  of  a  specific  test  and  associated  strength 
or  stiffness  criterion.  The  remainder  of  the  agencies  tend  to  use  grading  requirements. 


Elastic  Modulus 
21% 


Figure  2-6  Agencies  Method  of  Subgrade  Characterization 

Percentages  based  on  48  agencies  reporting  (updated  to  include  WSDOT). 


CHAPTERS 

PREDICTING  CRITICAL  STRAINS  AND  PAVEMENT  LIFE 

3.1  INTRODUCTION 

The  pavement  manager  is  concerned  with  the  life  remaining  in  any  given 
pavement  section.  A  powerful  tool,  then,  is  the  ability  to  estimate  the  pavement  life 
remaining  based  on  deflection  basin  parameters.  The  most  fundamental  of  these 
deflection  basin  parameters  are  the  Area  Parameter  and  Dq.  The  Area  Parameter  is  a 
normalized  width  of  the  deflection  basin  (refer  to  Figure  2-1)  and  Do  is  the  pavement 
surface  deflection  at  the  center  of  the  applied  load. 


P 


Figure  3-1  Typical  Deflection  Basin  (after  WSDOT  Pavement  Guide,  1995) 
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3.2  METHODOLOGY 

The  goal  of  the  research  discussed  in  this  chapter  is  to  develop  predictive 
relationships  that  can  predict  critical  strains  and  pavement  life  using  only  Dq,  Area 
Parameter  and  subrade  modulus  (Esg)-  With  such  relationships,  there  would  be  no  need 
to  know  the  layer  thicknesses,  AC  modulus,  or  base  modulus  to  estimate  the  critical 
strains  and  thereby  the  pavement  life  remaining.  These  relationships  would  also  be 
independent  of  the  pavement  temperature  because  they  would  be  based  on  data  generated 
at  a  wide  range  of  AC  stiffnesses  representing  a  wide  range  of  pavement  temperatures. 
3.2.1  Data  Generation 

In  order  to  generate  data  for  analysis,  a  matrix  of  different  pavement  scenarios 
was  established  (see  Table  3-1).  Three  parameters  were  chosen  to  be  variable  in  the 
pavement  scenarios  and  all  other  parameters  were  kept  constant.  The  variable  parameters 
were  subgrade  modulus,  asphalt  concrete  modulus  and  asphalt  concrete  thickness.  The 
constant  parameters  were: 


Poisson’s  Ratio,  AC  layer; 

0.35 

Poisson’s  Ratio,  Base  layer; 

0.40 

Poisson’s  Ratio,  Subgrade  layer: 

0.45 

Base  Modulus: 

207  MPa 

Base  Thickness: 

20  cm 

Load 

40  kN 

Tire  Pressure 

690  MPa 
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Table  3-1  shows  the  120  combinations  of  variable  parameters  and  indicates  which 
pavement  scenarios  (according  to  ‘case  number’)  correspond  to  which  combinations.  For 
example,  pavement  scenario  (case  number)  34  represents  a  Subgrade  Modulus  of  70 
MPa,  an  AC  Modulus  of  5520  MPa,  and  an  AC  Thickness  of  10  cm.  Initially  only  80 
pavement  scenarios  were  chosen.  Case  Numbers  81-120  were  added  because  additional 
data  points  were  needed  to  further  investigate  trends  that  were  identified  using  the  initial 
80  data  points. 

Each  pavement  scenario  was  analyzed  using  WSDOT’s  EVERSTRS  Program  to 
estimate  the  deflections  and  the  critical  strains  needed  to  predict  the  number  of  loads  to 
failure  for  both  the  rutting  criteria  and  the  AC  fatigue  criteria.  Figure  3-2  provides  a 
view  of  the  key  locations  of  concern  for  pavement  response.  Locations  2  and  4  are  the 
critical  strains  evaluated  in  this  thesis.  A  sample  EVERSTRS  output  can  be  found  in 
Appendix  A  along  with  a  compilation  of  the  critical  outputs  for  every  pavement  scenario. 
Note:  throughout  this  thesis,  the  negative  sign  normally  associated  with  vertical 
compressive  strain  in  the  subgrade  is  omitted.  A  negative  sign  indicates  a  case  in  which 
the  strain  at  the  top  of  the  subgrade  is  predicted  to  be  in  tension  instead  of  in 
compression. 

The  estimations  of  Dq,  Dj,  D2,  D3,  St  (horizontal  tensile  strain  at  bottom  of  AC 
layer),  and  (vertical  compressive  strain  at  the  top  of  the  subgrade  layer)  were 
transferred  to  a  spreadsheet  where  the  Area  Parameter,  Nf  for  Rutting,  and  Nf  for  AC 
Fatigue  could  be  calculated.  Different  parameters  were  plotted  against  each  other  to  see 
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if  any  relationships  could  be  seen.  The  results  were  very  promising  and  showed  a 
significant  relationship  between  key  measured  or  calculated  parameters  and  those  critical 
strains  required  for  computing  pavement  life.  Appendix  B  contains  the  graphs  that 
include  all  120  pavement  scenarios. 

The  pavement  scenarios  with  an  AC  thickness  of  40  cm  create  outlying  points  on 
some  of  the  graphs  presented  in  Appendix  B.  Looking  at  Figures  B-3  to  B-6  and  Figures 
B-8  to  B-1 1,  there  are  several  points  on  each  graph  that  have  a  combination  of  low  strain 
and  low  Area  Parameter  or  Dq.  These  points  that  do  not  fall  on  the  curve  are  the  very 
thick  40  cm  pavements. 

The  equations  used  to  calculate  loads  to  failure,  Nf  are: 

Rutting  Criteria  : 


where  Nf  =  allowable  number  of  80  kN  single  axle  loads  to  ensure  that  rutting  at  the 
pavement  surface  does  not  exceed  12.7  mm. 

8v  =  vertical  compressive  strain  at  the  top  of  the  subgrade 

AC  Fatigue  Criteria: 

Nj-  =  N^(iab)*  [shiftfactor ) 

where  log  uat)  =  14. 82  -  3. 291  -  0-  854 

where  Nf  =  number  of  80  kN  4  layer  (psi) 
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The  shift  factor  for  the  AC  fatigue  criteria  equation  is  generally  accepted  to  be  between  4 
and  10  (WSDOT  Pavement  Guide,  1995). 

3.2.2  Regression  Analysis 

Regression  creates  an  equation  which  predicts  the  variation  of  one  parameter 
based  on  the  values  of  one  or  more  other  parameters.  Appendix  C  provides  more 
detailed  information  on  regression  analysis.  The  equation  developed  through  linear 
regression  typically  follows  the  form  of: 

_y  =  + fljXj  +  a  2X2 + . 

where  y  is  the  dependent  variable  (the  predicted  value) 

Xj.  Xn  are  the  independent  variables  (the  known  values) 

ao.. 811  are  constants  determined  by  the  relationship  between  the  independent  and 
dependent  variables. 

The  relationship  is  easiest  to  visualize  when  there  is  one  independent  variable.  The 
equation  is  providing  the  “best  fit”  line  between  the  dependent  and  independent  variable. 
The  first  constant  is  the  y-intercept,  while  the  second  constant  is  the  slope  of  the  line. 
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Tire  with  specified 
load  and  pressure 


Legend 

1.  pavement  surface  deflection 

2.  horizontal  tensile  strain  at  bottom  of  AC 

3.  vertical  compressive  strain  at  top  of  base 

4.  vertical  compressive  strain  at  top  of  subgrad 


Figure  3-2  Pavement  Response  Locations  Used  in  Evaluating  Load  Effects 
(after  WSDOT  Pavement  Guide,  1995) 


Table  3-1  Case  Numbers  for  EVERSTRS  Calculations 
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There  are  several  important  descriptors  of  a  regression  equation  which  give  an 
indication  of  how  well  the  equation  will  predict  the  dependent  variable.  The  Coefficient 
of  Determination  is  annotated  as  R-Squared  or  R^.  R^  can  be  described  as  the  fraction  of 
the  variation  in  the  dependent  variable  that  is  explained  by  the  regression  equation  (Ryan, 
et.al.,  1985).  R^  can  range  from  0%  to  100%  with  the  latter  indicating  a  perfect  fitted 
line.  Another  important  descriptor  is  the  Root  Mean  Square  Error  (RMSE)  which  is  a 
measure  of  how  much  the  observed  dependent  variable  values  differ  from  the  average 
value  on  the  “best  fit”  line.  RMSE  is  the  standard  deviation  of  the  “best-fit”  or 
regression  equation  line  (Mahoney,  1994).  An  often  overlooked  descriptor  for  a 
regression  equation  is  the  number  of  observed  values  used,  designated  as  N.  Intuitively, 
the  larger  N  is,  the  more  weight  a  predictive  relationship  should  be  given. 

Regression  analysis  provides  an  equation  which  describes  a  straight  line  fit. 

Often,  variables  relate  to  each  other  in  a  curvilinear  manner  which  does  not  work  well 
with  a  direct  correlation  of  the  variables.  In  this  case,  a  transformation  is  done  on  one  or 
more  of  the  variables  so  that  they  relate  to  each  other  in  a  linear  relationship.  The  most 
common  transformations  are  the  log,  inverse,  square,  or  square  root  of  the  variable. 

All  of  the  data  generated  for  the  120  pavement  scenarios  was  entered  into  the 
statistical  software  program.  Minitab.  Regression  equations  were  computed  using 
different  transformations  of  the  variables  in  different  combinations.  The  equations  with 
the  highest  R^  and  lowest  RMSE  are  presented  in  the  following  section  while  the 
remainder  of  the  equations  developed  can  be  found  in  Appendix  D. 
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3.3  RESULTS 

3.3.1  Tensile  Strain  Restriction 

Note  that  the  equations  presented  in  the  next  two  sections  for  tensile  strain  at  the 
bottom  of  the  AC  layer  only  have  an  N  value  of  41 .  That  means  only  41  of  the  120 
pavement  scenarios  were  used  to  generate  these  equations.  This  is  because  the 

relationships  involving  Sj  were  discretely  curvilinear  only  for  Area  Parameters  greater 

than  500  (41  cases).  This  also  happens  to  be  the  region  in  which  the  tensile  strain  would 
tend  to  be  more  critical  than  the  vertical  compressive  strain  in  the  subgrade.  Figure  3-3 
shows  the  relationship  of  strain  at  the  bottom  of  the  AC  layer  to  Area  Parameter  for  a 
subgrade  modulus  of  35  MPa.  As  can  be  seen,  a  strain  of  200  would  predict  the  Area 
Parameter  to  be  approximately  440  or  620  depending  on  which  side  of  the  crest  one 
looks.  The  equations  discussed  above  used  only  those  points  to  the  right  of  the  crest, 
with  the  Area  Parameter  greater  than  or  equal  to  500.  Refer  to  Table  2-3  to  see  that  these 
are  relatively  strong  pavements.  These  are  the  pavements  for  which  AC  fatigue  is 
expected  to  be  more  critical  that  rutting  failure. 


Subgrade  Modulus  =  35  MPa 
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E  400 
%  300 
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0.00  200.00  400.00  600.00  800.00  1000.00 

Area  Parameter 


Figure  3-3  Dlustration  of  Area  Parameter  Limit  for  Prediction 
of  Strain  at  the  Bottom  of  the  AC  Layer 


3.3.2  Using  Area  Parameter 

The  best  equations  for  predicting  the  critical  strains,  8t  and  Sy,  using  Area 

Parameter,  A,  and  subgrade  modulus.  Ego,  as  calculated  on  Minitab  are: 

logSj  =  5.91-0. 00353A-0.8771ogEsQ 

r2  =  39  30/^  jyVjSE  =  0.1358  N  =  41 

and 

log  8y  =  7. 04  -  0. 003 3  9  ^  - 1. 34  log 

r2  =  9440/p  BMSE  =  0.1148  N  =  120 

where  A  is  in  mm  and  Eg©  is  in  MPa  and  is  calculated  using  one  of  the  closed  form 
equations  presented  in  Table  2-4.  These  equations  both  have  strong  R  values  showing 
strong  relationship  between  the  dependent  variable  and  the  independent  variables. 
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3.3.3  Using  Dn 

The  best  equations  for  predicting  the  critical  strains,  £t  and  £v,  using  deflection  at 

the  center  of  the  load.  Do,  and  subgrade  modulus.  Ego,  as  calculated  on  Minitab  are: 

log  8t  =  -4. 23  + 1. 80  log  Dq  +  0. 898  log 

=  97.6%  RMSE  =  0.0647  N  =  41 

and 

8^  =  -5 1 1  + 1. 97Do  +  0. 788Ese 

r2  =  96  0%  RMSE  =  163.9  N  =  120 

where  Dq  is  in  micrometers  and  Ego  is  in  MPa  and  is  calculated  using  one  of  the  closed 
form  equations  presented  in  Table  2-4.  These  equations  show  that  Dq  is  an  even  better 
predictor  for  critical  strains  than  Area  Parameter. 

3.3.4  Illustration  of  Equation  Results  Using  Sampling  of  Case  Numbers 

Thirteen  pavement  scenarios  from  Table  3-1  were  chosen  to  calculate  strains 
using  the  equations  presented  in  this  Chapter.  These  strains  were  then  compared  to  the 
strains  that  were  predicted  by  EVERSTRS.  The  equations  require  knowledge  of  Dq, 

Area  Parameter,  and  Ego-  Dq  and  Area  Parameter  were  taken  from  the  deflections 
estimated  by  EVERSTRS  (See  Table  A-1  in  Appendix  A  for  a  complete  listing  of  these 
deflections).  Ego  was  calculated  from  the  deflections  using  two  of  the  closed  form 
equations  presented  in  Table  2-4.  The  AASHTO  equation  was  used  in  conjunction  with 
the  deflection  at  two  feet  while  the  WSDOT  three  layer  equation  was  used  in  conjunction 
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with  the  deflection  at  three  feet.  The  other  equations  can  also  be  used  to  calculate  Esc 
from  deflection  measurements.  The  subgrade  modulus  calculated  from  the  deflection  at 
three  feet  produced  a  smaller  error  in  predicting  Sj  while  the  subgrade  modulus  calculated 
from  D2  resulted  in  a  smaller  error  in  predicting 

Table  3-2  presents  the  results  of  the  comparisons  for  the  tensile  strain  at  the 
bottom  of  the  AC  layer  and  Table  3-3  presents  the  results  of  the  comparisons  for  the 
vertical  compressive  strain  at  the  top  of  the  subgrade.  Table  3-2  only  contains  the  six 
case  numbers  which  had  an  Area  Parameter  greater  than  500  mm.  The  equation 
presented  by  Gomez- Achecar  and  Thompson  (1986)  for  Asphalt  strain  (see  Table  2-5) 
was  also  calculated  for  these  six  case  numbers  to  see  how  well  it  would  predict  the  tensile 
strain;  the  average  error  was  66%. 

For  the  case  numbers  presented  in  Tables  3-2  and  3-3,  the  Area  Parameter 
equation  appears  to  be  the  better  predictor  for  Sy  and  St.  While  these  errors,  ranging  from 
-7  to  +30  percent,  are  by  no  means  negligible,  they  are  certainly  reasonable  estimators 
considering  no  information  concerning  the  pavement  structure  and  layer  properties  is 
required. 

3.4  CASE  STUDY  COMPARISON 

Volume  3  of  the  WSDOT  Pavement  Guide  includes  four  case  studies  depicting 
actual  pavement  rehabilitation  projects.  Case  Study  Number  One  will  be  used  to 
compare  EVERSTRS  predicted  strains  with  the  strains  predicted  by  the  regression 
equations  presented  herein.  The  pavement  characteristics  at  each  core  location  are 
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presented  in  Table  3-4.  The  thicknesses  were  established  from  as-built  drawings  and 
core  samples.  The  asphalt,  base,  and  subgrade  moduli  were  back-calculated  using 
EVERCALC  with  no  stiff  layer  considered.  Two  core  locations  (MP  208.50  and  MP 
209.00)  were  not  used  in  this  comparison  because  the  root  mean  square  errors  were  large 
enough  to  consider  the  data  unreliable.  The  E^c  is  the  actual  backcalculated  modulus  for 
the  AC  at  the  insitu  field  temperatures. 

Table  3-2  Comparison  of  8t  as  Predicted  by  EVERSTRS  and  the  Predictive 
Equations  for  Selected  Case  Numbers 


Case 

St 

Ego  calculated  with  D2 

No. 

EVERSTRS 

St 

%  difference 

St 

%  difference 

Equation 

Equation 

using  Do 

using  A 

2 

361 

532 

+47 

659 

+83 

7 

303 

244 

-19 

214 

-29 

16 

43 

75 

+74 

33 

-23 

38 

176 

182 

+3 

173 

-2 

40 

23 

41 

+78 

19 

-17 

58 

161 

156 

-3 

174 

+8 

1 

Average  Error 

+30 

Average  Error 

+3 

Case 

St 

Ego  calculated  with  D3 

No. 

EVERSTRS 

St 

%  difference 

St 

%  difference 

Equation 

Equation 

using  Do 

using  A 

2 

361 

586 

+62 

600 

+66 

7 

303 

246 

-19 

212 

-30 

16 

43 

66 

+53 

38 

-12 

38 

176 

210 

+19 

151 

-14 

40 

23 

37 

+31 

20 

-13 

58 

161 

197 

+22 

139 

-14 

Average 

+28 

Average 

-3 

Error 

Error 

Table  3-3  Comparison  of  Ey  as  Predicted  by  EVERSTRS  and  the  Predictive 
Equations  for  Selected  Case  Numbers 
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Ev 

EVERSTRS 


Esg  calcxilated  with  D2 


%  difference 


Equation  using 

Do 


2081 


1091 


2288 


154 


1592 


911 


545 


-100 


964 


253 


837 


2178 


1636 


Average  Error 


Sv 

Equation  using 
A 


2106 


676 


2442 


85 


1988 


884 


415 


35 


1141 


303 


753 


2198 


2051 


Average  Error 


%  difference 


Sy 

EVERSTRS 


Esg  calculated  with  D3 


Ev 

Equation  using 
Do 

%  difference 

Ev 

Equation  using 
A 

2577 

2084 

-19 

1826 

733 

1091 

+49 

669 

2189 

2291 

+5 

2071 

119 

147 

+24 

103 

1745 

1604 

-8 

1490 

922 

922 

0 

680 

546 

554 

+1 

336 

58 

-110 

-290 

40 

1069 

992 

-7 

827 

371 

281 

-24 

214 

709 

891 

+26 

558 

2076 

2181 

+5 

1871 

1885 

1647 

-13 

1548 

%  difference 


Average  Error 


Average  Error 
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3.4.1  Prediction  of  Critical  Strains 

The  information  presented  in  Table  3-4  was  used  to  estimate  the  critical  strains 
using  EVERS TRS.  The  regression  equations  were  then  calculated  for  each  core  location. 
The  values  of  Dq  presented  in  the  WSDOT  manual  were  used  in  the  regression  equations. 
Area  Parameter  was  calculated  using  the  equation  presented  in  Table  2-2  using  the 
deflections  presented  in  the  WSDOT  manual.  The  deflection  values  were  already 
normalized  to  a  standard  load  of  40  kN  which  is  consistent  with  the  parameters  of  the 
predictive  equations.  If  these  values  had  not  already  been  normalized  for  load,  the 
procedure  developed  in  Chapter  4  would  have  been  used  to  normalize  Dq  and  Area 
Parameter.  The  subgrade  modulus  required  as  an  input  to  the  predictive  equations  was 
calculated  using  the  closed  form  equation:  Ego  =  -530  +  0.00877(P/D3) 
where  Ego  is  in  psi,  P  is  the  load  in  pounds,  and  D3  is  in  inches.  Conversions  were  made 
from  the  metric  units  used  in  this  thesis  to  the  U.  S.  Customary  units  for  this  equation. 
Once  Ego  was  calculated,  it  was  converted  back  to  MPa. 

Tables  3-5  and  3-6  present  the  results  for  the  strain  at  the  bottom  of  the  AC  layer 

(£,).  Tables  3-7  and  3-8  present  the  results  for  the  strain  at  the  top  of  the  subgrade  layer 

(Sy).  The  bold  values  in  Tables  3-5  and  3-6  indicate  that  the  equation  should  not  be  used 
for  this  core  location  because  the  Area  Parameter  is  less  than  500  mm.  Recall,  the 


predictive  equations  for  8t  are  valid  only  for  Area  Parameter  values  greater  than  500  mm 
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(see  Figure  3-4).  These  bold  values  are  not  included  in  the  calculation  of  the  average 
numerical  and  percent  differences  presented  at  the  bottom  of  the  tables. 

The  results  indicate  that  the  Dq  equation  better  predicts  8,  (-6%  error)  while  the 

Area  Parameter  equation  better  predicts  Sy  (+17%  error).  Again,  the  margin  of  error  is 

not  insignificant,  but  it  does  show  promise  as  a  pavement  management  tool. 

3.4.2  Prediction  of  Pavement  Life  Remaining 

The  critical  strains  predicted  in  the  last  section  (St  and  Sy)  are  amoung  the  inputs 
required  to  calculate  the  pavement  life  remaining  (measured  in  number  of  ESAL  loads). 
Refer  to  Section  3.2. 1  for  the  equations  to  calculate  pavement  life  remaining  (Nf)  for 
rutting  and  fatigue  failure  modes.  The  pavement  life  remaining  was  calculated  for 
WSDOT  Case  Study  Number  One  using  the  critical  strains  predicted  by  EVERSTRS  and 
the  critical  strains  predicted  by  the  predictive  relationships  presented  in  this  Chapter. 

For  each  mode  of  failure,  the  pavement  life  remaining  was  calculated  based  on 
the  critical  strains  predicted  by  EVERSTRS  and  the  predictive  relationship  (using  Dq  or 
A)  that  most  closely  matched  the  EVERSTRS  strain.  The  results  of  this  comparison 
were  very  promising  and  are  presented  in  Table  3-9  for  fatigue  failure  and  in  Table  3-10 
for  rutting  failure.  In  most  cases,  the  predicted  life  remaining  using  the  critical  strain 
predicted  by  the  predictive  equation  is  of  the  same  order  of  magnitude  as  the  predicted 
life  remaining  using  the  critical  strain  predicted  by  EVERSTRS. 


Table  3-5  Comparison  of  Results  for  Predicting  St  Using  Area  Parameter  &  E. 


Core  Using  Area  Parameter  and  Esq 

Location 

c  o  numerical  percent 

difference  difference 


EVERSTRS  I  EQUATION 


207.85 

287 

159 

-128 

-45 

208.00 

283 

170 

-113 

-40 

209.05 

480 

370 

-110 

-23 

209.40 

707 

379 

-328 

-46 

209.80 

261 

156 

-105 

-40 

210.00 

248 

207 

-41 

-17 

210.50 

137 

113 

-24 

-17 

211.00 

83 

53 

-30 

-36 

211.50 

81 

55 

-26 

-32 

212.00 

127 

65 

-62 

-49 

212.50 

94 

63 

-31 

-33 

Table  3-6  Comparison  of  Results  for  Predicting  8t  Using  Do& 


Core 

Location 

St 

EVERSTRS 

Using  Do . 

St 

EQUATION 

andEsG 

numerical 

difference 

percent 

difference 

207.85 

287 

246 

-41 

-14 

208.00 

283 

-42 

-15 

209.05 

480 

+64 

+13 

209.40 

707 

579 

-128 

-18 

209.80 

261 

234 

-27 

-10 

210.00 

248 

265 

+17 

+7 

210.50 

137 

181 

+32 

211.00 

83 

75 

-8 

-10 

211.50 

81 

74 

-7 

-9 

212.00 

127 

103 

-24 

-19 

212.50 

94 

91 

-3 

-3 

Average 

-22 

-6 

Table  3-8  Comparison  of  Results  for  Predicting  Ey  Using  Dq  &  E^g 


Core 

Location 

8v 

EVERSTRS 

Using  Do  i 

Sv 

EQUATION 

indEsG 

numerical 

difference 

percent 

difference 

207.85 

311 

492 

+181 

+58 

208.00 

256 

418 

+162 

+63 

209.05 

639 

1303 

+664 

+104 

209.40 

648 

1510 

+862 

+133 

209.80 

227 

387 

+160 

+70 

210.00 

402 

617 

+215 

+53 

210.50 

236 

265 

+29 

■B 

211.00 

0 

-108 

211.50 

51 

-10 

-61 

-120 

212.00 

98 

166 

■i"68 

+69 

212.50 

77 

97 

+20 

+26 

Average 

+199 

+33 
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Table  3-9  Comparison  of  Pavement  Life  Remaining  for  Fatigue  Failure  Criteria 


Core  Location 

NfLab 

No  Shift  Factor 
based  on 
EVERSTRS 
prediction 

NfLab 

No  Shift  Factor 
based  on 

predictive  equation 
prediction 

Percent 

difference 

207.85 

25,469 

42,300 

+66 

208.00 

32,569 

55,261 

+70 

209.05 

4,584 

3,037 

-34 

209.40 

3,095 

5,971 

+93 

209.80 

39,657 

56,807 

+43 

210.00  24,919  20,034  -20 

210.50  77,724  31,080  ^ 

211.00  911,725  1,272,692  +40 

211.50  1,338,522  1,802,222  +35 

212.00  454,203  904,942  +^ 

212.50  581,972  647,531  Tu 


Average  Difference 


+31 
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Table  3-10  Comparison  of  Pavement  Life  Remaining  for  Rutting  Failure  Criteria 


Core  Location 

Nf 

based  on 
EVERSTRS 
prediction 

Nf 

based  on 

predictive  equation 
prediction 

Percent 

difference 

207.85 

7,143,874 

8,916,733 

+25 

208.00 

17,098,133 

10,240,930 

o 

1 

209.05 

282,821 

92,428 

-67 

209.40 

265,628 

58,228 

-78 

209.80 

29,315,170 

16,511,743 

-44 

210.00 

2,259,895 

1,837,177 

-19 

210.50 

24,624,639 

78,955,256 

+221 

211.00 

819,854,957 

2,669,914,071 

+226 

211.50 

23,711,674,966 

4,748,106,725 

-80 

212.00 

1,267,547,032 

322,278,446 

-75 

212.50 

3,737,909,781 

492,475,622 

-87 

Average  Difference 

-2 

3.5  CONCLUSION 


The  regression  equations  developed  in  this  chapter  are  designed  to  provide  an 
approximate  estimate  of  the  critical  strains  in  the  pavement  structure.  These  critical 
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strains  can  then  be  used  to  predict  the  pavement’s  remaining  life.  These  estimates  are  not 
for  use  as  a  design  tool. 

Area  Parameter  and  Dq  must  be  normalized  to  a  40  kN  load  in  order  to  use  these 
equations  and  are  not  adjusted  for  temperature  as  presently  configured.  The  ranges  of 
AC  stiffness  used  in  the  data  generation  provides  for  a  range  of  temperatures.  The 

regression  equations  that  predict  St  can  only  be  used  with  pavement  sections  with  an 

Area  Parameter  greater  than  500  mm.  Pavements  with  an  Area  Parameter  less  than  500 
mm  often  fail  due  to  rutting,  although  fatigue  failure  may  also  control  for  these 
pavements.  Rutting  failure  is  not  typical  for  pavements  with  an  Area  Parameter  greater 


than  500  mm. 


CHAPTER  4 


LOAD  ADJUSTMENT 


4.1  INTRODUCTION 

When  testing  a  section  of  pavement,  two  very  important  factors  can  affect  the 
deflections  that  will  be  measured:  load  and  temperature.  As  the  load  increases,  the 
deflections  will  increase  for  the  same  pavement.  Deflections  also  increase  as  AC 
temperature  increases.  For  these  reasons,  test  results  must  be  adjusted  to  some  standard 
load  and  temperature  to  make  any  comparisons  valid.  This  chapter  describes  the  process 
used  to  develop  load  adjustment  factors  for  Area  Parameter  and  Dq.  The  process  used  to 
develop  a  temperature  adjustment  is  detailed  in  Appendix  F  because  temperature 
adjustment  is  not  required  with  the  equations  developed  in  this  thesis. 

The  load  adjustment  must  be  used  with  the  predictive  equations  developed  in 
Chapters  3  and  5.  The  temperature  adjustment  should  not  be  used  with  these  equations 
because  of  the  range  of  AC  stiffnesses  (hence  temperatures)  inherent  in  the  data  used  to 
develop  the  equations. 

A  load  adjustment  is  required  if  the  testing  will  use  other  than  a  standard  40  kN 
load.  Deflections  and  the  Area  Parameter  calculated  by  those  deflections  will  increase  as 
the  load  is  increased  for  the  same  pavement. 
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4.2  METHODOLOGY 
4.2.1  Data  Generation 

The  initial  120  pavement  scenarios  that  were  created  to  establish  relationships 
between  pavement  response  and  critical  strains  were  created  at  a  single  wheel  load  of  40 
kN.  In  order  to  develop  relationships  between  pavement  response  and  load,  10  pavement 
scenarios  were  chosen  to  run  additional  EVERSTRS  calculations.  The  pavement 
scenarios  were  chosen  over  a  range  of  pavement  thicknesses,  pavement  moduli,  and 
subgrade  moduli  and  are  listed  in  Table  4-1.  A  complete  listing  of  the  EVERSTRS 
output  values  for  these  10  pavement  scenarios  at  loads  of  20,  30,  40,  50,  60,  70,  and  80 
kN  can  be  found  in  Appendix  A. 

When  the  values  of  Area  Parameter  and  Dq  were  plotted  against  load,  as  can  be 
seen  in  Figures  4-1  and  4-2,  the  variation  of  Area  Parameter  and  Dq  were  approximately 
linear  with  changes  in  load.  Equations  relating  Area  Parameter  and  Do  with  load  were 
developed  using  the  regression  analysis  software  program  Minitab  and  are  discussed  in 


the  following  section. 


Table  4-1  Pavement  Scenarios  Used  to  Generate  Load  Adjustment  Factors 
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4.2.2  Regression  Analysis 

The  variation  for  both  Area  Parameter  and  Dq  was  estimated  as  a  function  of  load 
only.  In  mathematical  terms,  the  variation  of  Area  Parameter  and  Dq  due  to  load  is  the 
slope  of  the  “best  fit”  line  for  each  pavement  scenario.  The  slope  for  each  of  the 
pavement  scenarios  is  fairly  constant  for  Area  Parameter,  while  the  slope  of  the  line  for 
Do  can  be  seen  to  increase  with  increasing  y-intercept  values.  The  Area  Parameter 
variation  was  the  most  straightforward  to  determine. 

4.2.2. 1  Area  Parameter 

The  values  of  load  vs  Area  Parameter  for  each  pavement  scenario  was  analyzed 
using  the  Minitab  computer  program.  For  each  case,  a  regression  equation  was  generated 
for  Area  Parameter  with  Load  as  the  single  predictor.  The  equations  were  of  the  form: 

A  =  (intercept) +(slope)(Load) 

The  resulting  regression  equations  all  had  an  R  greater  that  90%  and  can  be  found  in 
Appendix  D.  The  intercept  ranged  in  value  from  201  mm  to  761  mm,  but  for  the 
purpose  of  establishing  a  load  adjustment  factor  the  intercept  is  not  important.  The  slope 
ranged  in  value  from  0.715  to  1.37.  The  differences  in  the  slope  seemed  to  be  random, 
so  an  average  slope  was  calculated  and  used  as  the  load  adjustment  factor  for  Area 
Parameter.  The  average  slope,  and  thus  the  load  adjustment  factor  is  1.069.  The 
equation  to  adjust  an  Area  Parameter  for  load  is: 


^adjusted  to  40kN  ^Load  (40  -  Load)( 1. 069) 
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As  an  example,  if  an  FWD  test  conducted  using  a  20  kN  load  produced  an  Area 
Parameter  of  450,  the  Area  Parameter  adjusted  to  40  kN  would  be  calculated  as  follows: 

A^j„3ted  to  401CN  =  450  +  (40  -  20)( 1. 069) 

A40.n  =  471.38 

If  the  actual  load  is  greater  than  40  kN,  Area  Parameter  will  be  shifted  to  a  lower  value, 
while  the  opposite  is  true  for  loads  less  than  40  kN.  Refer  to  Figure  4-1  for  a  visual 
representation  of  this  trend. 

4.2.2.2  Deflection  at  the  Center  of  the  Load 

Again,  the  values  of  Load  vs  Dq  were  analyzed  using  the  Minitab  computer 
program  to  find  the  regression  equation  describing  the  “best-fit”  line  for  each  pavement 
scenario.  These  equations  were  of  the  form: 

Dq  =  (intercept)  +  (slope)(Load) 

2 

The  regression  equations  had  excellent  R  values;  all  exceeded  98%  and  three  were 
perfect  fits  with  an  R  of  100%.  These  equations  can  be  found  in  Appendix  D.  The 
intercept  values  ranged  from  8  to  316  microns.  The  slope  values  ranged  from  3.76  to 
3 1 .0.  The  slope  values  increased  as  the  intercept  values  increased.  This  relationship 
created  the  opportunity  to  approximate  the  slope  (load  adjustment  factor)  rather  than  use 
an  average  as  was  done  for  the  Area  Parameter  load  adjustment  factor. 

The  slope  values  were  entered  into  Minitab  and  correlated  with  the  values  of  Do 
at  each  load.  Seven  regression  equations  were  developed  of  the  form: 


slope  =  (pintercept)  +  (pslope)(Do)L<„d 
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Pintercept  and  pslope  are  the  intercept  and  slope  from  the  equations  relating  the  slope 
pictured  in  Figure  4-2  to  Dq  at  each  load  The  goal  was  to  predict  slope  using  load  only, 
pintercept  and  pslope  were  correlated  with  Load  which  resulted  in  the  following 
regression  equations; 

— —  =  8.51  +  0. 97(Load)  for  all  loads 
pslope 

R^  =  99.9%  RMSE  =  0.5917  N  =  7 

which  can  be  simplified  as: 

~  (8.51  + 0.970( Load)) 
and 

pintercept  =  -1.27  +  0. 00758(Load)  for  loads  30-80  kN 

r2  =  97  4o/„  rMSE  =  0.0261 1  N  =  6 

pintercept  = -0.91  for  20  kN  load 

The  pintercept  equation  was  limited  to  the  30-80  kN  loads  because  the  value  for  the  20 
kN  load  was  severely  out  of  line  with  the  other  values.  The  regression  equation  for 
pintercept  which  included  the  20  kN  load  had  an  of  75.4%. 

An  example  to  see  how  these  equations  come  together  is  appropriate.  Given  a 
FWD  test  done  at  65  kN  and  a  measured  Dq  of  890  pm,  what  should  the  Dq  load 
adjustment  factor  be?  Given  the  load,  first  calculate  the  pintercept  and  pslope: 


pintercept  =  -1.27  +  0.00758(65  kN) 
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pintercept  =  -0.78 

~  (8.51  +  0.970(65  kN)) 
pslope  =  0.0140 

The  calculated  pintercept  and  pslope  are  then  inserted  into  the  equation  for  the  Dq  load 
adjustment  factor: 

slope  =  pintercept  +  pslope(Do)Load 
slope  =  -0.78 +  0.0140(890) 

Do  load  adjustment  factor  =  slope  =  1 1 .68 
The  adjusted  Dq  could  then  be  calculated  using  the  following  equation: 

(Po)40kN  “  (Po)Load  (40  -  Load)(Do)sF 

where  (Do)sf  is  the  Dq  load  adjustment  factor,  or  “shift  factor”.  The  adjusted  Dq  would 
be: 

(Do)40kN  =  (890)  + (40-65X11.68) 

(Do)40kN  =  598  pm 

The  procedure  to  adjust  Dq  for  load  is  more  complicated  that  that  for  Area 
Parameter.  In  summary  the  steps  are: 

1 .  Calculate  pintercept  and  pslope.  Input  required:  Load  used  in  testing. 

2.  Calculate  Dq  load  adjustment  factor.  Input  required:  Dq  measured  in  testing. 


3 .  Calculate  adjusted  Do . 
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4.3  RESULTS 

The  previous  section  described  the  process  to  adjust  Area  Parameter  and  Dq  to  a 
standard  load  of  40  kN.  How  well  these  parameters  can  be  adjusted  significantly  impacts 
their  use  in  predicting  critical  strains  necessary  to  calculate  life  left  in  the  pavement 
structure.  Table  4-2  shows  how  well  Area  Parameter  was  shifted  for  load  for  the  10 
pavement  scenarios  used  to  generate  the  initial  data.  Column  3  shows  Area  Parameter 
calculated  by  EVERSTRS  for  that  Case  Number  at  40  kN  load,  column  5  shows  Area 
Parameter  calculated  by  EVERSTRS  for  that  Case  Number  at  the  load  indicated  in 
column  2,  and  column  4  shows  Area  Parameter  adjusted  to  40  kN  using  the  process 
outlined  in  this  chapter.  On  average,  the  shifted  value  of  Area  Parameter  is  within  1%  of 
the  actual  Area  Parameter  measured  at  40  kN.  Table  4-3  shows  how  well  Dq  was  shifted 
for  load  for  the  same  10  pavement  scenarios.  On  average,  the  error  in  shifting  Dq  to  40 
kN  is  6.4%.  These  comparisons  use  only  theoretical  data  created  by  the  EVERSTRS 
computer  program;  however,  results  should  be  comparable  for  FWD  data  since 
EVERSTRS  was  created  based  on  FWD  data. 


Table  4-2  Error  in  Adjusting  Area  Parameter  to  40  kN 


Table  4-2  Error  in  Adjusting  Area  Parameter  to  40  kN,  continued 


Case  No. 


Percent 

Difference 


605.7 

605.7 


605.71 


605.71 


264.36 


26 


.36 

.36 


601.33 

605.981 


602.99 


598.84 


593.75 


588.00 


260.27 


262.86 


265.22 


265 

265 


579.95 


595.29 


613.68 


620.22 


625.82 


630.76 


238.89 


2.17| 

5.91 


87.07| 

97.95 


352.01 


352.01 


460.24 


460.24 


460.24 


460.24 


460.24 


460.24 


25 


250.86 


250.86 


250.86 


250.86 


250.86 


347.00 


460.52 


458.81 


248.59 


252.52 


253.57 


254.34 


254.87 


379.07 


446.64 


472.16 


482.85 


492.59 


501.57 


223.91 


237.90 


263.21 


74.9 

86.4 


297.63 


0.72% 


0.04% 


0.45% 


1.97% 


2.92% 


1.55% 


0.57% 


0.33% 


0.50% 


65 

80 

308.63 

0.57% 

70 

20 

352.01 

327.84 

0.79% 

70 

30 

352.01 

351.75 

341.06 

0.07% 

70 

50 

352.01 

351.03 

361.72 

0.28% 

70 

60 

352.01 

349.28 

370.66 

0.78% 

2.1 

1.8 


.63% 


0.27% 


0.27% 


0.06% 


0.31% 


1.08% 


1.39% 


1.60% 


61.53% 


1.03% 


Total  Error 


Average  Error 
(excluding  40  kN  load) 


Table  4-3  Error  in  Adjusting  Dq  to  40  kN 


Case 

Niunber 


Do 

calculated 
at  Load 


906.10 


1259.69 


1589.13 


1901.36 


2200.60 


2489.10 


2768.73 


544.45 


794.55 


1037.40 


1274.13 


1505.75 


1732.64 


1955.37 


697.14 


009.56 

307.67 


1594.13 


1871.01 

2139.40 


2400.48 


121.01 


158.93 


196.54 


233.95 


271.16 


308.22 


266.05 


378.26 


484.11 


584.75 


681.05 


11 

62 


Do 

calculated 
at  40  kN 


1589.13 


1589.13 

1589.13 


1589.13 


1589.13 


1589.13 


1589.13 


1037.40 


1037.40 


1037.40 


1037.40 


1037.40 


1037.40 


1037.40 


1307.67 


1307.67 
1307.67' 
1307.67' 


1307.67 

1307.67' 


1307.67 

158.93 


158.93 


158.93 


158.93 


158.93 


158.93 


158.93 


484.11 


484.11 


484.11 


484.11 


484.11 


Do 

adjusted 
to  40  kN 


1537.20 


1584.20 


1589.13 


1576.76 


1557.15 


1534.01 


1509.14 


916.40 


995.38 


1037.40 


1059.55 


1070.62 

1074.56 


1073.60 


1178.50 


1267.56 


1307.67 

1323.42 


1326.38 

1321.61 


1311.95 
23.69 


142.76 


158.93 


170.98 


180.11 


186.88 


438.50 


468.41 


484.11 


491.09 


493.17 


491.99 


488.43 


Percent 

Difference 


3.27% 


0.31% 


0.00% 


0.78% 


5.03% 


11.66% 


2.13% 


3.20% 


3.58% 


3.49% 


9.88% 


3.07% 


0.00% 


1.20% 


1.43% 


1.07% 


0.33% 


22.17% 


10.17% 


0.00% 


7.58% 


13.33% 


17.59% 


20.55% 


9.42% 


3.24% 


0.00% 


1.44% 


1.87% 


50 

70 

50 

80 

Table  4-3  Error  in  Adjusting  Dq  to  40  kN,  continued 


Case 

Number 

Load 

(kN) 

Do 

calculated 
at  Load 

Do 

calculated 

at40kN 

Do 

adjusted 
to  40  kN 

Percent 

Difference 

55 

40 

239.58 

239.58 

239.58 

0.00% 

55 

50 

294.01 

251.35 

4.91% 

55 

60 

347.57 

259.67 

8.39% 

55 

70 

400.33 

239.58 

265.33 

10.75% 

55 

80 

452.38 

268.78 

12.19% 

65 

20 

365.45 

574.91 

609.13 

5.95% 

65 

30 

574.91 

596.52 

3.76% 

657.40 


730.38 

795.94 


855.65 


200.13 


279.93 


353.73 


422.68 


607.47 


845.31 


1199.06 


1528.74 


1840.40 


2138.11 

2424.191 


2700.61 


431.79 


545.54 


636.36 


712.45 


779.64 


839.36 


893.75 


574.91 


574.91 


57 


574.91 


353.73 


353.73 


353.73 


353.73 


353.73 


353.73 


353.73 


1528.74 


1528.74 


1528.74 


1528.74 


1528.74 

1528.74' 


1528.74 


636.36 


636.36 


636.36 


636.36 


636.36 


636.36 


636.36 


527.71 


505.62 


343.93 


353 

357 


357.74 


355.68 


351.83 


1432.85 

1507.45 


1528.74 


1526.49 


1513.40 


1494.59 


1472.66 


723.01 


680.17 


636.36 


201 


531.99 


505.13 


8.21% 


12.05% 


15.69% 


8.03% 


2.77% 


1.13% 


0.55% 


0.54% 


0.00% 


0.15% 


3.67% 


13.62% 


6.88% 


11.65% 


16.40% 


20.62% 


383.68% 


6.39% 


Total  Percentage  Error 


Average  Percentage  Error 
(excluding  40  kN  loads) 


65 


4.4  CONCLUSION 

The  load  shift  factors  presented  in  this  chapter  were  developed  using  the  data 
created  for  this  thesis.  Any  shift  factor  can  be  used  with  the  various  equations  presented 
in  this  thesis  as  long  as  the  pavement  engineer  has  confidence  in  the  system  being  used. 
The  critical  point  is  to  adjust  pavement  responses  to  a  standard  load  (40  kN)  before  using 
the  relationships  presented  throughout  this  thesis.  Temperature  shift  factors  are 
presented  in  Appendix  F;  however,  these  shift  factors  should  not  be  used  with  the 
equations  presented  in  this  thesis  because  temperature  variation  was  included  in  the  data 
used  to  develop  the  equations. 


CHAPTER  5 


PREDICTING  OVERLAY  REQUIREMENTS 

5.1  INTRODUCTION 

Chapter  3  outlined  the  steps  taken  to  estimate  life  remaining  in  a  pavement  using 
measurable  pavement  responses  and  parameters  that  can  be  directly  calculated  from  those 
responses,  namely;  deflection  at  the  center  of  the  load  (Dq),  Area  Parameter  (A),  and 
subgrade  elastic  modulus  (Esc)  Another  powerful  tool  is  the  ability  to  predict  the 
overlay  required  for  a  pavement  system  based  on  those  parameters  and  the  number  of 
equivalent  single  axle  loads  (ESALS)  the  pavement  needs  to  be  designed  for.  This 
chapter  discusses  the  research  conducted  to  develop  relationships  that  can  predict  the 
overlay  required  based  on  these  readily  available  parameters. 

5.2  METHODOLOGY 
5.2.1  Data  Generation 

A  layered  elastic  analysis  was  performed  for  each  of  the  120  pavement  scenarios 
originally  shown  in  Table  3-1.  The  analysis  reported  in  this  chapter  was  done  using 
EVERPAVE  4.0.  The  input  parameters  that  were  held  constant  for  this  analysis  were: 
Design  Tire  Load  40000  N 

Tire  Pressure  690  kPa 

Dual  Spacing  0  cm  (i.e.,  a  single  tire  was  used) 


Fatigue  Shift  Factor,  New  AC 


10 
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Fatigue  Shift  Factor,  Old  AC  10 

Seasonal  Variation  All  seasons  at  25  degrees  Celsius 

Lane  Distribution  F actor  1 . 0 

Minimum  Overlay  0.5  cm  (a  value  of  0  is  not  allowed  by  the 

program) 

Each  case  was  analyzed  for  total  design  traffic  (80  kN  ESALS)  of  1,  5,  and  10  million. 
The  overlay  thickness  required  and  the  damage  levels  for  fatigue  damage  on  the  new  AC, 
fatigue  damage  of  the  old  AC  and  rutting  damage  on  the  subgrade  were  recorded.  These 
values  can  be  found  in  Appendix  A,  Table  A-2.  The  overlays  were  significantly  greater 
than  typically  expected  and  the  reasons  for  this  will  be  described  in  Section  5.3. 

5.2.2  Regression  Analysis 

Pavement  scenarios  that  resulted  in  an  overlay  requirement  of  0.5  cm  only  due  to 
the  program  requirement  that  some  minimum  overlay  be  considered  were  removed  for 
the  regression  analysis  of  the  data.  This  was  done  because  these  values  tended  to  skew 
the  data.  Pavements  with  a  wide  range  of  deflection  responses  will  need  the  minimum 
overlay  dictated  by  EVERPAVE. 

The  results  from  the  EVERPAVE  analyses  were  entered  into  the  Minitab 
computer  program  to  run  the  regression  analysis.  Regression  analysis  was  done  using 
Do,  Area  Parameter,  Ego,  and  ESALS  as  independent  variables.  The  overlay  thickness 
required  was  the  dependent  variable.  All  of  the  equations  calculated  using  these 
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variables  in  various  transformations  can  be  found  in  Appendix  D.  The  best  equation 
found  was; 

Overlay  =  -541  +  84.91ogDo+  82.41ogA  +  54.31ogEgG+  14.51ogESALS 
R^  =  95.7%  RMSE  =  2.239  N  =  327 

where  overlay  is  in  cm,  Dq  and  Area  Parameter  are  in  pm  and  mm  respectively.  Ego  is  in 
MPa  and  ESALS  is  in  millions. 

5.3  RESULTS 

When  the  overlay  equation  was  compared  with  the  overlay  predicted  by 
EVERPAVE  for  the  327  cases  analyzed,  the  overall  results  were  very  promising.  On 
average  there  was  a  13%  difference  from  the  overlay  predicted  from  the  equation 
compared  to  the  overlay  predicted  by  EVERPAVE  (refer  to  Appendix  A,  Table  A-3  for 
the  complete  listing  of  the  comparisons).  The  absolute  difference  on  average  was  1.4 
cm.  This  shows  that,  given  constant  parameters,  overlay  can  be  predicted  fairly  well 
from  Do,  Area  Parameter,  Ego,  and  ESALS  without  knowing  the  composition  of  the 
pavement  system.  Once  the  equation  was  compared  with  some  of  the  case  studies  in  the 
WSDOT  Pavement  Guide,  however,  it  became  clear  that  some  of  the  simplifying 
assumptions  created  an  overprediction  of  the  overlay.  Two  Case  Studies  discussed  in  the 
WSDOT  Pavement  Guide  were  used  to  compare  overlay  results. 

5.3.1  Case  Study  Number  One 

Case  Study  Number  One  in  the  WSDOT  Pavement  Guide  covers  the  section  of 
State  Route  395  from  Mile  Post  207.81  to  Mile  Post  212.67.  The  pavement  guide 
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determined  the  required  overlay  using  three  methods.  The  first  method  is  WSPMS 
SCOPER  which  is  a  component  analysis  approach  partially  based  on  the  Asphalt 
Institutes  Component  Analysis  procedure.  The  second  method  is  DARWin,  the 
computerized  version  of  the  pavement  design  models  in  AASHTO’s  “Guide  for  Design 
of  Pavement  Structures  1993”.  The  third  and  final  method  is  using  EVERPAVE  4.0,  the 
mechanistic-empirical  overlay  design  procedure  developed  by  WSDOT. 

Table  5-1  shows  the  comparison  of  the  results  from  the  three  methods  used  in  the 
pavement  guide  along  with  the  results  using  the  overlay  equation  presented  in  this 
chapter.  No  temperature  adjustment  was  made  for  this  comparison.  The  last  four  core 
locations  may  have  shown  an  overlay  requirement  due  only  to  the  minimum  requirement 
imposed  by  the  computer  program.  For  this  reason,  these  core  locations  were  not 
included  in  the  computation  of  the  average  difference  between  EVERPAVE  and  the 
predictive  equation.  The  average  overlay  thickness  difference  between  the  EVERPAVE 
results  and  the  Equation  result  is  20.6  cm.  The  two  primary  factors  which  may  be 
contributing  to  this  difference  are  the  seasonal  effects  and  the  load  configuration  chosen. 

The  overlay  equation  is  based  on  EVERPAVE  data  which  was  executed  with  no 
seasonal  variation  and  an  average  temperature  of  25  degrees  Celsius.  This  would 
represent  a  climate  about  the  equivalent  of  Miami,  Florida.  In  order  to  quantify  the 
effect  on  the  results,  the  first  nine  core  locations  were  recalculated  with  EVERPAVE;  the 
only  change  being  the  seasonal  effects  were  held  to  the  constant  25  degrees  Celsius.  The 
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Table  5-2  Overlay  Discrepancy  Based  on  Seasonal  Variation 


Core 

Location 

EVERPAVE 
w/  seasonal  variation 
(cm) 

Recalculation  @  25  "C 
EVERPAVE 
(cm) 

Difference 

(cm) 

207.85 

5.0 

11.5 

6.5 

208.00 

4.0 

11.5 

7.5 

208.50 

1.0 

6.0 

5.0 

209.00 

9.5 

16.0 

6.5 

209.05 

10.0 

16.5 

6.5 

209.40 

11.0 

19.0 

8.0 

209.80 

3.5 

10.0 

6.5 

210.00 

6.0 

11.0 

5.0 

210.50 

3.5 

10.0 

6.5 

Another  significant  difference  results  from  the  way  the  tire  load  was  modeled.  In 
developing  the  data  used  to  develop  the  overlay  equation,  the  load  was  modeled  as  a 
single  40  kN  load  rather  than  dual  tires  loaded  to  20  kN  each  and  spaced  35.6  cm  apart 
which  is  more  typical  of  the  assumption  used  in  the  original  WSDOT  case  study.  The 
single  load  was  used  because  EVERPAVE  analyzes  pavement  much  faster  with  a  single 
load  rather  than  a  dual  load.  Several  of  the  original  case  scenarios  were  recalculated  with 
a  dual  load,  all  other  parameters  were  held  constant.  The  results  are  shown  in  Table  5-3. 


Table  5-3  Overlay  Discrepancy  Based  on  Load  Configuration 


Case  No 

Single  Tire  Overlay 
(cm) 

Dual  Tire  Overlay 
(cm) 

Difference 

(cm) 

1 

34.0 

19.0 

15.0 

10 

26.0 

11.5 

14.5 

31 

12.5 

0 

>12.5 

74 

17.0 

3.5 

13.5 

81 

35.0 

20.5 

14.5 

85 

33.5 

19.0 

14.5 

94 

28.0 

14.0 

i 

14.0 

100 

25.0 

11.0 

14.0 

108 

26.0 

11.5 

14.5 

120 

22.0 

9.0 

13.0 

The  average  difference  due  to  the  load  configuration  is  14.2  cm.  This  added  with  the 
difference  due  to  temperature  (6.4  cm)  combines  for  a  total  average  difference  of  20.6 
cm.  This  is  also  the  average  difference  that  was  found  between  the  EVERPAVE  results 


and  the  overlay  equation  results  in  Table  5-1. 
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5.3.2  Case  Study  Number  Two 

Case  Study  Number  Two  in  the  WSDOT  Pavement  Guide  covers  the  section  of 
Washington  State  Route  7  from  Mile  Post  0.00  to  Mile  Post  16.82.  Again,  the  pavement 
guide  determined  the  overlay  required  using  three  methods. 

Table  5-4  shows  the  comparison  of  the  results  from  the  three  methods  used  in  the 
pavement  guide  along  with  the  results  using  the  overlay  equation  presented  in  this 
chapter.  No  temperature  adjustment  was  made  for  this  comparison. 

The  average  difference  between  the  EVERPAVE  results  and  the  overlay  equation 
results  is  24.06  cm.  This  is  a  3.5  cm  increase  over  the  average  difference  for  Case  Study 
Number  One.  Again,  the  two  primary  factors  contributing  to  this  difference  are  the 
seasonal  effects  and  the  tire  spacing  chosen. 

5.4  CONCLUSION 

The  overlay  equation  developed  from  the  360  cases  calculated  by  EVERPAVE 
shows  a  high  correlation  of  the  data.  This  indicates  that  the  overlay  required  could  be 
reasonably  estimated  using  only  Do,  Area  Parameter,  Esc,  and  ESALS.  This  equation 
would  be  an  extremely  useful  tool  for  prioritizing  pavements  for  rehabilitation.  It  would 
also  give  an  early  estimation  of  the  extent  of  overlay  required.  The  particular  equation 
developed  in  this  thesis,  however,  severely  overestimates  the  overlay  required  when 
compared  with  previous  case  studies  presented  in  the  WSDOT  Pavement  Guide.  The 
two  factors  contributing  to  this  overestimation  are  the  seasonal  temperature  and  load 
configuration  used  when  developing  the  data.  A  similar  overlay  equation  developed 
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using  standard  seasonal  variation  and  load  configuration  for  a  region  would  provide  more 

reasonable  results.  The  consistency  of  the  overlay  equation  results  developed  in  this 

chapter  are  very  promising,  indicating  further  development  would  be  worthwhile. 

Table  5-4  Comparison  of  Overlay  Results  -  Case  Study  Number  Two 
(No  Temperature  Adjustment) 


Core  WSPMS 

Location  overlay  (cm) 


0.23  4.9 


AASHTO 
overlay  (cm) 


0 


EVERPAVE 
overlay  (cm) 


0.5 


Equation 
overlay  (cm) 


16.5 


Table  5-4  Comparison  of  Overlay  Results  -  Case  Study  Number  Two 
(No  Temperature  Adjustment),  continued. . . 
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Core 

Location 

WSPMS 
overlay  (cm) 

AASHTO 
overlay  (cm) 

EVERPAVE 
overlay  (cm) 

Equation 
overlay  (cm) 

8.78 

10.1 

11.4 

4.5 

28.5 

9.48 

10.7 

13.0 

4.0 

34.6 

9.98 

13.4 

14.7 

8.0 

35.0 

10.58 

11.3 

13.7 

6.0 

36.2 

10.98 

10.7 

10.7 

7.0 

29.2 

11.63 

10.7 

10.7 

7.5 

30.7 

12.18 

9.8 

14.7 

5.5 

36.8 

12.53 

9.4 

11.4 

5.0 

31.5 

■1 

9.8 

7.4 

1.5 

26.6 

13.53 

9.4 

8.4 

1.0 

23.0 

14.53 

7.3 

7.6 

0.5 

21.7 

15.08 

7.6 

9.1 

3.5 

24.8 

15.68 

8.2 

15.0 

8.0 

30.2 

16.03 

0 

5.1 

1.0 

23.3 

16.55 

0 

3.0 

0.5 

18.8 

CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  CONCLUSIONS 

General  conclusions  concerning  the  research  presented  in  the  previous  chapters 
are  included  here. 

•  Critical  strains  (Sj  and  Sy)  can  be  estimated  using  only  Esg  and  Area  Parameter  or  Dq. 

From  these  critical  strains,  pavement  life  remaining  can  be  estimated. 

•  Overlay  required  can  be  estimated  using  Area  Parameter,  Dq,  Ego,  and  ESALS  given 
consistent  seasonal  variation  and  load  configuration. 

•  Area  Parameter  and  Do  must  be  adjusted  to  a  standard  load  (40  kN)  to  use  the 
equations  presented  in  this  thesis.  Temperature  adjustments  are  not  required. 

•  The  equations  presented  are  estimations  and  should  be  used  at  the  pavement 
management  level,  not  the  design  level.  More  rigorous  methods  are  required  for 
design  purposes. 

•  The  equations  developed  and  presented  in  this  thesis  are  summarized  in  Table  6-1 . 
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Table  6-1  Summary  of  Equations  Developed  in  this  Thesis 


Equation 

RMSE 

N 

log  St  =  5.91  -  0.00353A  -  0.8771ogEsG 

89.3 

0.1358 

41 

log  St  =  -4.23  +  l.SlogDo  +  0.8981ogEsG 

97.6 

0.0647 

41 

log  Sv=  7.04  -  0.00339A  -  1.341ogEsG 

94.4 

0.1148 

120 

Sv  =  -511  +  1.97Do  +  0.788Esg 

96.0 

163.9 

120 

overlay  =  -541  +  84.91ogDo  +  82.41ogA  +  54.31ogEsG  + 

95.7 

2.239 

327 

14.51ogESALS 

6.2  RECOMMENDATIONS 

The  pavement  scenarios  used  to  create  the  data  in  this  thesis  ranged  in  AC 
thickness  from  2.5  cm  to  40  cm.  The  very  thick  40  cm  pavements  do  not  correlate  as 
well  as  the  pavements  20  cm  and  thinner.  The  predictive  equations  could  be 
reformulated  without  the  very  thick  pavements  which  should  increase  the  R  and 
decrease  the  RMSE.  Of  course  this  would  also  limit  the  range  of  pavements  the 
equations  could  be  applied  to. 

Additional  case  study  comparisons  should  be  made  for  both  the  critical  strain  and 
the  overlay  predictive  equations.  A  wide  range  of  pavements  should  be  compared  with 
these  equations;  however,  the  pavements  should  be  within  the  range  of  variable  pavement 
parameters  used  in  this  research  (see  Table  3-1). 
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The  critical  strain  equations  could  be  recalculated  using  an  array  of  pavement 
scenarios  that  include  variations  in  base  moduli  and  thickness.  Another  variation  that 
could  be  attempted  would  be  to  limit  the  AC  moduli  to  the  equivalent  for  25 '  C.  This 
would  create  a  regression  equation  that  would  need  to  be  adjusted  for  temperature.  A 
comparison  could  be  made  between  these  types  of  equations  to  determine  which  is  a 
better  predictor  of  critical  strains. 

Further  research  is  required  for  the  overlay  equation.  An  overlay  equation  could 
be  developed  using  a  more  standard  seasonal  variation  and  a  dual  load  configuration. 
Several  seasonal  variations  could  be  tried  to  see  if  there  is  a  correlation  between  average 
annual  air  temperature  and  additional  overlay  required. 

The  equations  presented  in  this  thesis  represent  many  hours  of  data  generation 
and  manipulation.  As  they  exist  now,  they  can  be  used  as  an  approximate  estimation  of 
the  pavement  condition  and  the  extent  of  overlay  required.  It  is  hoped  that  further 
research  will  result  in  better  predictive  relationships  that  can  be  used  with  confidence  for 
prioritizing  pavement  rehabilitation  projects. 
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APPENDIX  A 
DATA 

Figures  A-1  and  A-2  show  a  sample  output  for  the  computer  program 
EVERSTRS.  Each  of  the  120  pavement  scenarios  was  analyzed  using  EVERSTRS.  The 
key  results  are  listed  in  Table  A-1.  This  includes  the  deflections  at  each  sensor,  the  Area 
Parameter  calculated  from  these  deflections  and  the  strains  at  the  bottom  of  the  asphalt- 
concrete  layer  and  at  the  top  of  the  subgrade.  Table  A-2  presents  the  EVERPAVE 
results  for  the  120  pavement  scenarios  at  three  different  levels  of  ESALs  required.  The 
overlay  required  and  the  percent  of  damage  accumulated  for  rutting  and  AC  fatigue  are 
shown  in  this  table.  Table  A-3  is  a  comparison  of  the  EVERPAVE  prediction  of  overlay 
required  (as  shown  in  Table  A-2)  and  the  overlay  predicted  by  the  overlay  equation 


presented  in  Chapter  5. 


Layered  Elastic  Analysis  by  Everstrs  4,0 

Title:  Input! 

No  of  Layers:  3 

No  of  Loads:  I 

No  of  X-Y  Evaluation  Points:  4 

Layer 

Poisson's 

Thickness 

Moduli(l) 

Ratio 

(cm) 

(MPa) 

1 

J5 

10.000 

690.00 

2 

.40 

20.000 

207.00 

3 

.45 

35.00 

Load  No 

X-Position 

Y-Position 

Load  Pressure 

Radius 

(cm) 

(cm) 

(N) 

(kPa) 

(cm) 

1 

.00 

.00 

40000.0  690.00 

13.584 

Location  No:  I 

X-Position  (cm):  .000 

Y-Position  (cm): 

.000 

Normal  Stresses 

Z-Position 

Layer 

Sxx 

Syy 

Szz  Syz 

Sxz 

Sxy 

(cm) 

(kPa) 

(fcPa) 

(kPa)  (kPa) 

(kPa) 

(kPa) 

9.999 

1 

410.93 

410.93 

-403.32  .00 

.00 

.00 

20.000 

2 

68.01 

68.01 

-172.81  .00 

.00 

.00 

30.001 

3 

-5.18 

-5.18 

-66.44  .00 

.00 

.00 

.001 

1 

-1246.43 

-1246.43 

-690.00  .00 

.00 

.00 

Normal  Strains  tad  Deflettions 

Z-Position 

Layer 

Exx 

Eyy 

Ezz  Ux 

Uy 

Uz 

(cm) 

(10^-6) 

(!(r-6) 

(10^-6)  (microns) 

(microns) 

(microns) 

9.999 

1 

591.69 

591.69 

1001.41  .00 

.00 

1258.04 

20.000 

2 

531.07 

531.07 

1097.69  .00 

.00 

1120.12 

30.001 

3 

772.77 

772.77 

1764.99  .00 

.00 

1015.54 

.001 

I 

-824.18 

-824.18 

264.50  .00 

.00 

1302.61 

Principal  Stresses  and  Strains 

Z-Posilion 

Layer 

SI 

S2 

S3  El 

E2 

E3 

(cm) 

(W*a) 

(kPa) 

(kPa)  (10^-6) 

(i(r-6) 

(I(r-6) 

9.999 

I 

^3J2 

410.93 

410.93  -1001.41 

591.69 

591.69 

20.000 

2 

-172.81 

68.01 

68.01  -1097.69 

531.07 

531.07 

30.001 

3 

-66.44 

-5.18 

-5.18  -1764.99 

772.77 

772.77 

.001 

1 

-1246.43 

-1246.43 

-690.00  -824.18 

-824.18 

264.50 

Location  No:  2 

X-Position  (cm):  .000 

Y-Positioa  (cm): 

30.500 

Normal  Stresses 

Z-Position 

Layer 

Sxx 

Syy 

Szz  Syz 

Sxz 

Sxy 

(cm) 

(kPa) 

(kPa) 

(kPa)  (kPa) 

(kPa) 

(kPa) 

9.999 

1 

15.31 

-120.64 

-15.67  -61.75 

.00 

.00 

20.000 

2 

36.58 

-3.41 

-26.21  -5824 

.00 

.00 

30.001 

3 

-4J2 

-14.02 

-2921  -15.01 

.00 

.00 

.001 

1 

-206.56 

-13.19 

.00  -.01 

.00 

.00 

Normal  Slratns  and  Dcflecdons 

Z-Position 

Layer 

Exx 

Eyy 

Ezz  Ux 

Uy 

Uz 

(cm) 

(10^-6) 

(lir-6) 

(10^-6)  (microns) 

(microns) 

(microns) 

9.999 

I 

91J3 

-174.66 

30.71  .00 

27.86 

810.70 

20.000 

2 

233.94 

-36.51 

-190.70  .00 

7U5 

802.28 

30.001 

3 

432.44 

30.54 

-598.82  .00 

131.89 

773.49 

.001 

1 

•292.67 

85.65 

111.47  .00 

-8926 

803.22 

Principal  Stresses  and  Strains 

E3 

Z-Position 

Layer 

SI 

S2 

S3  El 

E2 

(cm) 

(kPa) 

(kPa) 

(kPa)  (10^-6) 

(10^-6) 

(10^-6) 

9.999 

1 

-149.20 

12.89 

15  J1  -230.53 

86.59 

91.33 

20.000 

2 

-74.15 

36.58 

44.53  -514.95 

233.94 

287.73 

30.001 

3 

-38.43 

-4.80 

-4J2  -980.90 

412.62 

432.44 

.001 

1 

-206.56 

-13.19 

.00  -292.67 

85.65 

m.47 

Location  No:  3 

X-Position  (an):  .000 

Y-Positiem  (cm): 

61.000 

Figure  A-1  Sample  EVERSTRS  Output  Sheet,  page  1 


Layered  Elastic  Analysis  by  Everstrs  4.0 


Normal  Stresses 

Z-Position 

Layer 

Sxx 

Syy 

Szz 

Syz 

Sxz 

Sxy 

(cm) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

9.999 

1 

-7.75 

-26.67 

-2.57 

-14.55 

.00 

.00 

20.000 

2 

9.34 

-13.54 

-6.92 

-14.07 

.00 

.00 

30.001 

3 

-2.80 

-10.08 

-9.45 

-6.87 

.00 

.00 

.001 

1 

-60.84 

32.73 

.00 

.00 

.00 

.00 

Normal  Strains  and  Deflections 

Z-Position 

Layer 

Exx 

Eyy 

Ezz 

Ux 

Uy 

Uz 

(cm) 

(l{r-6) 

(10'-6) 

(I0-.6) 

(microns) 

(microns) 

(microns) 

9.999 

1 

3.61 

-33.43 

13.73 

.00 

2.20 

517.79 

20.000 

2 

84.65 

-70.06 

-25.33 

.00 

51.64 

517.20 

30.001 

3 

171.15 

-130.50 

-104.46 

.00 

104.40 

513.27 

.001 

1 

-104.78 

78.30 

14.26 

.00 

-63.92 

516.30 

Principal  Stresses  and  Strains 

E3 

Z-Position 

Layer 

SI 

S2 

S3 

El 

E2 

(cm) 

(kPa) 

(kPa) 

(kPa) 

(10^-6) 

(10^-6) 

(icr-6) 

9.999 

1 

-33.52 

-7.75 

4.27 

^6.81 

3.61 

27.12 

20.000 

2 

-24.68 

9.34 

-145.46 

50.07 

84.65 

30.001 

3 

-16.65 

-2.89 

-2.80 

^2.49 

167.54 

171.15 

.001 

I 

-60.84 

.00 

32.73 

-104.78 

14.26 

78.30 

Location  No;  4 

X-Position  (cm):  .000 

Y-Position  (cm): 

91.500 

Normal  Stresses 

Z-Position 

Layer 

Sxx 

Syy 

Srz 

Syz 

Sxz 

Sxy 

(cm) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

(IcPa) 

9.999 

1 

-5J5 

-7J8 

-.84 

-3.13 

.00 

.00 

20.000 

2 

1.99 

-11.14 

•229 

-3.64 

.00 

.00 

30.001 

3 

-1.42 

-5.70 

-3.23 

-3.05 

.00 

.00 

.001 

1 

-20.78 

32.95 

.00 

.00 

.00 

.00 

Normal  Strains  and  Deflections 

2-Position 

Layer 

Exx 

Eyy 

Ezz 

Ux 

Uy 

Uz 

(cm) 

(10^-6) 

(10^-6) 

(10^-6) 

(microns) 

(microns) 

(microns) 

9.999 

I 

-3.45 

-7.60 

5.19 

.00 

-3.15 

354 J5 

20.000 

2 

35.54 

•5323 

6.64 

.00 

32.52 

354.85 

30.001 

3 

74.25 

-103.04 

-.74 

.00 

67.94 

355.72 

.001 

1 

-46.83 

58  JO 

-6.17 

.00 

-42.85 

354J7 

Principal  Stresses  aid  Strains 

E3 

Z-Position 

Layer 

SI 

S2 

S3 

El 

E2 

(cm) 

(kP«) 

(kPa) 

(kPa) 

(10^-6) 

(l(r-6) 

ii(r-6) 

9.999 

1 

-8.63 

•525 

.42 

-10.06 

-3.45 

7.65 

20.000 

2 

-12.44 

-.98 

1.99 

-62.06 

15.47 

35.54 

30.001 

3 

-7.75 

-1.42 

-1.17 

-188.24 

74.25 

84.46 

.001 

1 

-20.78 

.00 

32.95 

-46.83 

-6.17 

58.30 

Figure  A-2  Sample  EVERSTRS  Output  Sheet,  page  2 
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Table  A-1  Data  Output  from  EVERSTRS 


Case 

D1 

D2 

D3 

A 

Strain 

H.  Strain 

No. 

Subgrade  top 

AC  bottom 

(pm) 

(pm) 

(pm) 

(mm) 

(10-6) 

(10-6) 

1 

1589.13 

875.9 

527.88 

348.66 

455.09 

2576.61 

360.66 

2 

1302.61 

803.22 

516.3 

354.27 

502.60 

1764.99 

591.69 

3 

961.25 

661.62 

477.74 

353.15 

569.67 

942.49 

400.24 

4 

674.14 

463.52 

379.63 

313.3 

604.44 

393.84 

160.37 

5 

1499.90 

522.61 

349.67 

467.78 

2394.77 

471.07 

6 

1165.17 

778.10 

510.11 

354.68 

535.78 

1519.22 

522.98 

7 

799.12 

611.82 

461.53 

350.2 

628.58 

733.26 

303.46 

8 

512.18 

400.41 

343.89 

293.88 

682.78 

279.26 

113.29 

9 

1408.41 

838.84 

518.03 

349.96 

483.91 

2189.20 

485.38 

10 

1037.40 

748.18 

504.83 

355.68 

572.80 

1258.73 

409.25 

11 

667.21 

552.20 

437.64 

343.61 

683.07 

541.52 

207.72 

12 

397.60 

338.50 

302.41 

267.4 

746.22 

187.01 

72.53 

■a 

1307.67 

826.50 

515.24 

350.43 

505.98 

1942.84 

427.18 

■a 

911.64 

706.32 

497.48 

357.4 

614.63 

993.75 

291.18 

15 

553.08 

484.42 

403.38 

329.81 

732.54 

379.09 

131.69 

16 

311.56 

280.51 

258.56 

235.64 

795.04 

119.48 

43.41 

17 

1194.72 

810.49 

514.34 

351.87 

535.28 

1654.59 

335.04 

18 

786.87 

648.90 

482.64 

357.82 

660.01 

744.02 

192.44 

19 

453.36 

413.54 

359.99 

306.82 

775.59 

252.89 

78.97 

20 

245.14 

228.93 

216.09 

201.85 

831.21 

73.76 

24.86 

21 

1104.64 

479.38 

259.87 

167.68 

379.51 

1866.10 

408.53 

22 

906.84 

458.03 

260.61 

171.09 

422.70 

1289.17 

587.44 

23 

675.37 

176.24 

485.73 

696.44 

380.95 

24 

493.57 

288.08 

217.54 

168.78 

516.76 

296.10 

146.549 

25 

1041.85 

474.05 

259.02 

168.22 

391.47 

1745.03 

479.68 
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Table  A-1  Data  Output  from  EVERSTRS,  Cont. 


Case 

Do 

D1 

D2 

D3 

A 

Strain 

H.  Strain 

No. 

Subgrade  top 

AC  bottom 

(pm) 

(pm) 

(pm) 

(pm) 

(mm) 

(10-6) 

(10-6) 

26 

797.47 

451.25 

260.53 

171.9 

457.30 

1112.70 

502.18 

27 

543.03 

369.51 

251.84 

177.89 

551.09 

544.08 

283.17 

28 

357.66 

248.99 

200.87 

162.63 

605.07 

212.49 

102.92 

29 

972.94 

471.58 

258.27 

168,49 

407.44 

1598.27 

472.03 

30 

697.39 

440.02 

261.66 

173.32 

496.95 

921.73 

384.48 

31 

441.87 

335.98 

244.31 

178.79 

614.35 

403.52 

191.55 

32 

268.16 

210.83 

179.68 

151.94 

682.62 

144.08 

65.91 

33 

893.44 

470.98 

258.36 

168.8 

430.01 

1416.57 

404.87 

34 

601.96 

420.08 

262.94 

175.99 

542.80 

727.43 

269.41 

35 

359.51 

296.49 

230,32 

176.67 

673.93 

284.45 

120.72 

36 

205.18 

175.09 

155.88 

137.14 

745.93 

93.26 

39.6 

37 

804.04 

468.12 

260.09 

169.62 

460.60 

1203.06 

311.98 

38 

511.36 

389.44 

261.26 

179.6 

593.78 

545.75 

176.25 

39 

291.04 

254.49 

209.89 

169.47 

727.48" 

191.68 

72.3 

40 

158.93 

143.22 

131.83 

119.88 

794.85 

58.31 

22.8 

41 

801.3 

253.06 

125.27 

'  81.32 

311.78 

1239.62 

447.87 

42 

658.84 

255.42 

127.04 

82.23 

348.36 

865.12 

583.69 

43 

497.16 

235.22 

131.29 

85.37 

403.27 

474.4r 

365.21 

180.63 

122.07 

87.86 

429.49 

205.72 

135.12 

45 

752.01 

253.82 

125.14 

81.36 

322.49 

1165.35 

484.84 

46 

566.24 

257.24 

128.02 

82.4 

381.96 

749.63 

483.07 

47 

383.6 

222.81 

133,06 

87.26 

469.83 

373.25 

265.47 

48 

261.57 

155.86 

115.44 

87.31 

519.41 

149.99 

93.78 

49 

696.08 

256.77 

125.02 

81.33 

337.38 

1069.39 

457.86 

50 

484.11 

255.84 

130.34 

82.91 

421.64 

622.06 

361.7 
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Table  A-1  Data  Output  from  EVERSTRS,  Cont. 


Case 

Do 

D1 

D2 

D3 

A 

Strain 

H.  Strain 

No. 

Subgrade  top 

AC  bottom 

(urn) 

(pm) 

(pm) 

(pm) 

(mm) 

(10-6) 

(10-6) 

51 

301.99 

204.54 

132.82 

89.7 

538.17 

278.77 

177.1 

52 

187.63 

131.99 

105.45 

84.17 

606.48 

103.31 

59.88 

53 

630.38 

261.49 

125.27 

81.28 

359.06 

947.52 

383.09 

54 

408.75 

248.23 

133.89 

84.35 

468.79 

491.98 

249.46 

55 

239.58 

181.79 

128.8 

91.44 

605.71 

198.36 

110.77 

56 

138.9 

109.74 

93.12 

78.21 

683.36 

67.99 

36.04 

57 

557.3 

265.31 

126.68 

81.35 

389.03 

803.99 

290.09 

58 

340.41 

233.08 

136.94 

87.16 

522.73 

370.83 

161.48 

59 

190.44 

156.94 

120.42 

90.83 

669.00 

135.35 

66.16 

60 

105.11 

89.88 

79.88 

70.08 

746.28 

43.23 

20.84 

61 

618.17 

129.21 

59.79 

39.87 

255.42 

751.72 

477.39 

62 

508.4 

141.02 

60.17 

39.8 

529.61 

581.42 

63 

388.98 

141.88 

65.35 

40.56 

294.79 

354.14 

313.01 

116.48 

67.71 

44.39 

353.37 

130.27 

126.78 

574.91 

131.6 

59.66 

39.82 

264.36 

709.22 

487.87 

66 

425 

146.17 

60.76 

39.67 

315.03 

460.60 

467.92 

67 

286.48 

136.39 

68.37 

41.66 

392.42 

233.80 

251.72 

68 

202.85 

99.74 

65.83 

45.67 

435.50 

96.60 

86.52 

69 

525.83 

136.16 

59.44 

39.75 

277.30 

651.64 

445.59 

353.73 

149.32 

62.54 

39.58 

352.01 

383.16 

343.18 

■■ 

216.77 

126.51 

70.74 

43.63 

460.43 

176.03 

165.37 

72 

138.25 

84.14 

61.6 

45.65 

524.04 

67.67 

54.82 

468.5 

142.78 

59.28 

39.64 

296.75 

577.35 

364.81 

■■ 

290.96 

148.06 

65.85 

39.91 

397.39 

304.04 

233.06 

75 

166.5 

113.21 

71.03 

45.93 

531.72 

126.53 

102.46 

92 


76 

98.12 

69.8 

55.49 

43.86 

609.72 

45.34 

32.93 

77 

405.97 

149.48 

59.78 

39.44 

324.32 

489.95 

271.88 

78 

236.48 

141.29 

69.81 

41.29 

451.10 

230.50 

149.21 

79 

129.02 

98.2 

68.49 

47.73 

602.57 

87.48 

60.9 

80 

71.86 

57.09 

48.38 

40.44 

685.52 

29.38 

19.08 

■■ 

1748.75 

920.5 

531.73 

343.7 

435.47 

3115.38 

-178.28 

■■ 

1669.6 

896.58 

530.03 

346.47 

444.47 

2836.18 

141.74 

B 

1693.88 

897.84 

528.31 

345.78 

440.13 

2989.49 

31.68 

84 

525.39 

347.95 

452.33 

2684.78 

316.38 

85 

1643.31 

877.46 

523.76 

346.76 

444.46 

2872.60 

195.37 

86 

1527.39 

856.52 

520.62 

348.41 

461.98 

2522.43 

408.63 

1591.29 

862.57 

519.58 

346.89 

450.36 

2745.93 

294.28 

1446.9 

845.23 

517.11 

348.5 

476.09 

2325.13 

416.23 

89 

1528.74 

853.85 

516.74 

346.82 

460.24 

2582.61 

324.75 

90 

1351.85 

837.05 

515.43 

348.99 

496.69 

2076.35 

363.85 

91 

1202.4 

492.55 

258.88 

165.69 

363.88 

2234.75 

-51.34 

92 

1156.5 

485.3 

259.47 

166.74 

370.66 

2046.30 

220.54 

93 

1174.29 

485.58 

258.57 

166.44 

367.15 

2163.91 

100.23 

94 

1112.39 

478.69 

258.77 

167.37 

377.40 

1951.88 

347.03 

95 

1146.63 

479.12 

257.71 

166.87 

370.44 

2092.00 

221.90 

96 

1062.89 

474.3 

257.83 

167.65 

386.39 

1840.68 

408.84 

97 

1112.43 

474.8 

256.73 

166.98 

375.71 

2005.19 

295.14 

1001.11 

473.28 

257.19 

167.76 

400.34 

1696.95 

400.70 

473.75 

256.09 

166.96 

385.29 

1885.43 

312.35 

100 

474.5 

257.43 

167.99 

421.34 

1511.81 

342.67 
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Table  A-1  Data  Output  from  EVERSTRS,  Cont. 


Case 

Do 

D1 

D2 

D3 

A 

Strain 

H.  Strain 

No. 

Subgrade  top 

AC  bottom 

(pm) 

(pm) 

(pm) 

(pm) 

(mm) 

(10-6) 

101 

863.02 

253.68 

124.6 

80.91 

300.29 

1470.95 

53.84 

102 

836.13 

252.9 

124.94 

81.12 

304.92 

1354.16 

285.67 

103 

847.08 

252.24 

124.59 

81.02 

302.57 

1434.70 

156.59 

104 

804.76 

252.12 

124.84 

81.19 

310.55 

1299.79 

371.01 

105 

830.39 

250.96 

124.41 

81.04 

305.06 

1393.65 

242.31 

106 

766.91 

252.8 

124.63 

81.17 

318.54 

1229.49 

406.28 

107 

250.72 

124.14 

80.98 

309.55 

1338.54 

293.28 

108 

ihh 

255.98 

124.45 

81.09 

331.50 

1133.86 

384.48 

109 

766.46 

252.69 

123.92 

80.86 

318.25 

1258.33 

299.04 

110 

652.74 

261.41 

124.59 

81.02 

351.56 

1008.87 

321.93 

111 

59.87 

39.9 

247.39 

885.45 

131.99 

112 

59.83 

39.88 

250.38 

818.59 

334.29 

HU 

59.81 

39.88 

249.01 

867.80 

198.48 

617.2 

127.76 

59.74 

39.85 

254.84 

788.99 

388.35 

115 

636.36 

125.93 

59.72 

39.84 

250.86 

845.42 

256.82 

116 

584.85 

129.81 

59.6 

39.79 

261.48 

747.68 

402.98 

117 

615.98 

126.81 

59.59 

39.77 

254.47 

812.79 

290.51 

118 

540.72 

134.22 

39.71 

272.71 

370.45 

119 

581.91 

129.59 

59.42 

39.69 

261.80 

287.35 

120 

485.1 

141.08 

59.09 

39.61 

290.62 

613.07 

304.48 
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Table  A-2  EVERPAVE  Output 


E  SG  ESALS  Overlay 


Rutting 

Damage 


(|im) 


1589 


1303 


961 


674 


1500 


1165 


799 


512 


1408 


1037 


667 


1308 


912 


553 


195 

787' 


453 


105 

907 


675 


494 


1042 


797 


543 


(mm) 


455 


503 


570 


604 


468 


536 


MPa 


35 


35 


35 


35 


35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

35 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

70 

1 

34.00 

0.536 

0.275 

0.527 

0.289 

0.698 

0.266 

20.00 

0.921 

0.081 

33.00 

0.372 

0.478 

29.00 

0.271 

0.505 

21.00 

0.264 

0.480 

4.00 

0.000 

0.459 

31.00 

0.205 

0.710 

26.00 

0.079 

0.705 

16.00 

0.025 

0.700 

11.00 


26.50 


20.00 


5.00 


28.50 


26.00 


27.001 


23.50 


15.00 


27.50 


22.50 


20.00 


8.50 


0.000 


0.008 


0.000 


0.000 


0.981 


0.930 


0.689 


0.240 


0.090 


0.020 


0.061 


0.004 


0.000 


0.848 


0.700 


0.727 


0.970 


0.506 


0.465 


0.910 


0.970 


0.951 


0.211 


0.911 


0.961 


0.932 


0.930 


0.976 


0.968 


0.958 


0.988 


0.911 


0.895 


0.946 


0.869 


0.802 


0.971 

0.995 

0.910 

0.916 

0.922 

0.973 

0.933 

0.927 

0.920 

0.617 

0.929 

0.555 

0.941 

0.555 

0.952 

0.512 
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Case 

Do 

No. 

(pm) 

37 

804 

38 

511 

39 

291 

41 

801 

42 

659 

43 

497 

45 

752 

46 

566 

47 

384 

49 

696 

50 

484 

51 

302 

53 

630 

54 

409 

55 

240 

57 

557 

58 

340 

59 

190 

miQQim 

HBI 

63 

389 

65 

575 

^■1 

WtEM 

im^l 

69 

526 

70 

354 

71 

217 

73 

469 

74 

291 

75 

167 

A  E  SG  ESALS 
(mm)  MPa  (millions) 


461  70 


594  70 


727  70 


312  140 


140 


140 


140 


140 


140 


Overlay  Fatigue 
Damage 
(cm)  New  AC 


23.00 


27.00 


24.50 


22.00 


25.50 


0.005 


25.50 


20.5 

10.5 


24.0 

18.5 


7.00 


21.50 


15.50 


24.00 


21.50 


24.50 


24.50 


19.50 


9.00 


23.00 


17.00 


5.50 


0.991 


0.937 


0.680 


0.239 


0.000 


0.003 


0.000 


0.924 


0.954 


0.649 


0.219 


0.092 


0.010 


0.047 


0.003 


0.000 


Fatigue 
Damage 
Old  AC 

Rutting 

Damage 

0.964 

0.682 

0.945 

0.581 

0.967 

0.397 

0.497 

0.249 

0.226 

0.103 

0.960 

0.290 

0.970 

0.297 

0.987 

0.319 

0.991 

0.215 

0.985 

0.208 

0.976 

0.205 

0.988 

0.193 

0.927 

0.171 

0.915 

0.150 

0.962 

0.221 

0.945 

0.181 

0.167 


0.907 


0.921 


0.915 


0.998 


0.916 


0.950 


0.911 


0.034 


51 

55 


0.037 


0.036 


0.040 


0.035 


0.034 


0.029 


96 


Case 

Do 

A 

ESG 

No. 

(urn) 

(mm) 

MPa 

77 

406 

324 

280 

78 

236 

451 

280 

81 

1749 

435 

35 

82 

1670 

444 

35 

83 

440 

35 

84 

452 

35 

85 

444 

35 

86 

462 

35 

87 

■n^ 

450 

35 

88 

1447 

476 

35 

1529 

460 

35 

90 

1352 

497 

35 

91 

1202 

364 

70 

92 

1157 

371 

70 

93 

1174 

367 

70 

94 

1112 

377 

70 

95 

1147 

70 

96 

1063 

70 

97 

1112 

376 

70 

98 

1001 

400 

70 

99 

1064 

385 

70 

100 

925 

70 

101 

863 

[IIIIIQjQ 

140 

102 

836 

305 

140 

103 

847 

mm 

140 

104 

805 

HB 

140 

105 

830 

305 

140 

106 

767 

319 

140 

107 

806 

310 

140 

108 

716 

332 

140 

ESALS 

Overlay 

(millions) 

(cm) 

1 

20.50 

1 

14.00 

1 

35.00 

1 

34.50 

1 

34.50 

1 

33.50 

1 

33.50 

1 

32.50 

1 

32.00 

1 

30.50 

1 

30.00 

1 

1 

1 

29.50 

1 

29.00 

1 

28.00 

1 

30.00 

1 

28.50 

1 

29.50 

1 

27.50 

1 

27.50 

1 

25.00 

1 

29.00 

1 

28.00 

1 

27.00 

1 

26.50 

1 

28.00 

1 

1 

1 

26.00 

Fatigue 
Damage 
New  AC 


0.003 


0.000 


0.610 


0.569 


0.369 


0.264 


0.209 


Rutting 

Damage 


0.957 


0.895 


0.779 


0.2101 


0.116 


0.073 


0.953 


0.956 


0.763 


0.906 


0.971 


0.267 


0.274 


0.468 


0.497 


.715 

.685' 


0.900 


0.809 


0.894 


0.774 


0.421 


0.468 


0.939 


0.977 


0.986 


0.972 


0.977 


0.950 


0.041 


0.037 


0.955 


0.934 


0.960 


0.988 


0.982 


0.928 


0.993 


0.942 


0.999 


0.985 


0.748 


0.807 


0.633 


0.668 


0.529 


0.568 


0.561 


0.642 


0.218 


0.233 


0.460 

0.998 

0.211 

0.315 

0.949 

0.199 

0.189 

0.940 

0.161 

0.102 

0.937 

0.174 

97 


Table  A-2  EVERPAVE  Output,  Cent. 


Overlay 


(pm) 


766 


653 


660 


643 


649 


617 


636 


585 


616 


541 


582 


485 


A 

ESG 

(mm) 

MPa 

318 

140 

352 

280 

249 

280 

255 

280 

251 

280 

261 

280 

254 

280 

280 

280 

291 

280 

455 

35 

503 

35 

570 

35 

604 

35 

mm\ 

mm\ 

629 

35 

683 

35 

484 

35 

573 

35 

683 

35 

746 

35 

506 

35 

615 

35 

733 

35 

535 

35 

660 

mm 

776 

35 

Fatigue 
Damage 
Old  AC 


Rutting 

Damage 


26.00 

0.062 

0.929 

0.017 

0.929 

0.964 

0.433 

27.00 

0.911 

0.437 

26.00 

0.923 

25.00 

0.975 

27.00 

0.927 

25.50 

0.969 

26.50 

0.178 

0.938 

0.094 

0.944 

0.055 

0.924 

0.014 

0.949 

0.877 

0.429 

0.807 

0.433 

34.00 

0.928 

29.50 

0.972 

HlfQIgg] 

41.00 

0.744 

0.750 

0.440 

0.737 

18.50 

0.352 

40.00 

0.969 

35.00 

0.160 

0.962 

25.00 

0.078 

0.956 

5.00 

0.000 

0.945 

38.50 

0.985 

32.00 

0.022 

0.981 

21.00 

0.001 

0.943 

35.00 

0.027 

0.985 

28.50 

0.000 

0.931 

16.00 

0.000 

0.963 

0.178 


0.207 


0.042 


0.040 


0.055 


0.055 


0.037 


0.039 


0.032 


0.035 


0.037 


0.988 


0.942 


0.258 


0.988 


0.387 


0.852 


0.820 


0.684 


0.705 


0.622 


0.722 


0.659 


0.544 
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Case 

Do 

A 

No. 

(urn) 

(mm) 

21 

1105 

380 

22 

907 

423 

23 

486 

24 

517 

25 

391 

26 

797 

457 

27 

543 

551 

28 

358 

605 

29 

973 

407 

30 

■BBI 

31 

nog 

32 

268 

683 

33 

893 

430 

34 

602 

543 

35 

360 

674 

37 

804 

461 

38 

594 

41 

mmi 

659  348 


497  403 


566  382 


MPa 


70 


70 


70 


70| 


70 


Overlay 

(cm) 


39.00 


36.00 


32.50 


29.50 


37.00 


33.00 


25.00 


5 

26.50 

5 

14.00 

5 

37.00 

34.50 

mam 

31.50 

5 

29.00 

5 

35.00 

5 

31.00 

5 

22.50 

5 

5.50 

5 

35.50 

5 

30.50 

5 

20.50 

5 

0.50 

Fatigue 
Damage 
New  AC 


0.971 


0.993 


0.623 


0.338 


0.157 


0.098 


0.000 


0.109 


0.018 


0.001 


0.020 


0.000 


0,964 


0.963 


0.967 


0.286 


.314 

.153 


0.123 


0.000 


Fatigue 
Damage 
Old  AC 

Rutting 

Damage 

0.466 

0.518 

0.485 

0.548 

0.318 

0.360 

0.090 

0.080 

0.628 

0.663 

0.666 

0.118 

0.988 

0.422 

0.981 

0.418 

0.975 

0.414 

0.961 

eeees 

0.937 

0.961 

0.3351 

0.949 

0.327 

0.940 

0.244 

0.453 

0.162 

0.424 

0.153 

0.245 

0.091 

0.963 

0.209 

0.975 

0.230 

0.970 

0.238 

0.957 

0.137 

0.157|| 
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Table  A-2  EVERPAVE  Output,  Cont. 


101 


Case 

Do 

No. 

(urn) 

1 

1589 

2 

Bigg 

3 

HKSI 

4 

674 

5 

1500 

6 

1165 

7 

799 

8 

512 

9 

1408 

10 

11 

12 

398 

13 

1308 

mm 

mm 

553 

mm 

1195 

18 

787 

19 

IBS 

HEI 

Bi^ 

907 

23 

675 

24 

25 

1042 

26 

797 

moBk 

543 

358 

29 

973 

30 

697 

31 

442 

32 

268 

33 

893 

34 

602 

mi 

!BI 


(mm) 


455 


503 


570 


604 


468 


536 


629 


683 


484 


573 


683 


746 


5 


6 


733 


535 


660 


776 


380 


423 


486 


517 


391 


457 


551 


605 


407 


497 


614 


683 


MPa 


35 


Overlay 

(cm) 


47.00 


44.00 


39.00 


34.50 


45.50 


41.50 


34.00 


23.50 


45.50 


40.50 


30.50 


10.50 


.5 


26.00 


40.50 


33.00 


21.00 


44.00 


41.50 


37.50 


34.00 


42.00 


38.00 


30.00 


23.50 


43.00 


38.00 


28.00 


7.50 


41.50 


35.00 


23.50 


Fatigue 
Damage 
New  AC 


0.943 


0.924 


0.966 


0.941 


0.756 


0.583 


0.506 


0.139 


0.146 


0.033 


0.001 


0.000 


0.944 


0.834 


0.662 


0.724 


0.979 


0.362 


0.179 


0.110 


0.029 


0.135 


0.028 


0.003 


Fatigue 
Damage 
Old  AC 

Rutting 

Damage 

0.450 

0.857 

0.491 

0.937 

0.849 

BESBli 

0.263 

0.923 

0.902 

0.994 

0.848 

0.949 

0.408 

0.377 

0.971 

0.677 

0.964 

0.671 

0.958 

0.664 

0.948 

0.653 

0.945 

0.494 

0.966 

Bl^ 

0.979 

0.976 

0.525 

0.995 

0.553 

0.433 

bce^^be^ 

0.464 

0.430 

0.342 

0.324 

0.113 

0.088 

0.938 

0.565 

0.281 

0.125 

0.951 

0.329 

0.945 

0.327 

mm 

0.323 

0.980 

0.974 

0.275 

0.995 

0.256 
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Table  A-2  EVERPAVE  Output,  Cont. 
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Case 

Do 

A 

No. 

(pm) 

(mm) 

70 

354 

352 

71 

217 

460 

73 

297 

74 

397 

75 

167 

532 

77 

406 

324 

78 

236 

451 

82 

83 

1694 

84 

1599 

mm 

85 

1643 

86 

462 

87 

450 

88 

1447 

HQI 

89 

1529 

BIBI 

90 

1352 

497 

91 

1202 

364 

92 

1157 

371 

93 

1174 

367 

1112 

95 

1147 

96 

1063 

386 

97 

1112 

376 

98 

1001 

400 

385 

100 

■ifflHni 

101 

863 

■■i 

102 

836 

305 

Esg  ESALS 
(MPa)  (millions) 


280  10 


280  10 


280  10 


280  10 


Overlay 

Fatigue 

Fatigue 

Rutting 

Damage 

Damage 

Damage 

(cm) 

New  AC 

Old  AC 

33.00 

0.186 

0.984 

0.027 

23.00 

0.193 

0.978 

37.00 

0.108 

0.942 

0.022 

30.50 

0.023 

0.989 

0.023 

19.50 

0.007 

0.019 

33.50 

0.020 

HOB 

0.026 

27.00 

0.000 

1.000 

0.023 

14.50 

0.000 

0.946 

0.015 

48.50 

0.980 

0.401 

0.804 

47.50 

0.989 

47.00 

46.00 

0.848 

48.00 

0.598 

0.978 

0.680 

47.50 


43.50 


46.00 


45.00 


45.50 


44.00 


43.50 


42.50 


0.465 


0.298 


0.076 


0.966 


0.967 


0.569 


0.463 


0.287 


0.200 


0.124 


0.993 


0.948 


0.958 


0.396 


0.436 


0.851 


0.965 


0.958 


0.9851 


0.978 


0.987 


0.948 


0.991 


0.3991 


0.652 


0.436 


0.455 


0.385 


0.411 


0.331 


0.344 


0.133 
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Table  A-2  EVERPAVE  Output,  Cont. 


Case 

Do 

A 

Esg 

ESALS 

Overlay 

Fatigue 

Fatigue 

Rutting 

No. 

Damage 

Damage 

Damage 

(pm) 

(mm) 

(MPa) 

(millions) 

(cm) 

New  AC 

Old  AC 

103 

847 

303 

140 

10 

42.00 

0.939 

0.810 

0.159 

311 

140 

10 

40.00 

0.976 

0.988 

0.192 

105 

305 

140 

10 

43.00 

0.549 

mEEs 

0.117 

106 

767 

319 

140 

10 

41.50 

0.443 

0.121 

140 

10 

43.50 

mE^E 

0.948 

0.084 

108 

716 

140 

10 

41.00 

mES^ 

0.968 

0.094 

109 

766 

318 

140 

10 

41.50 

0.115 

0.966 

0.081 

110 

653 

140 

10 

0.056 

0.983 

0.100 

111 

660 

IflQQ 

280 

10 

0.963 

0.396 

0.028 

112 

643 

280 

10 

40.50 

0.435 

0.029 

113 

280 

10 

39.50 

0.871 

0.037 

114 

617 

255 

280 

10 

38.00 

0.972 

0.993 

0.040 

115 

636 

251 

mg 

10 

40.50 

0.558 

0.027 

116 

585 

261 

mo^ 

10 

39.50 

0.421 

0.026 

117 

616 

254 

280 

10 

ieqss 

0.257 

0.962 

0.020 

118 

541 

273 

280 

10 

mm] 

0.174 

0.995 

0.023 

119 

582 

262 

280 

10 

39.00 

0.106 

0.970 

0.020 

120 

485 

291 

280 

10 

36.00 

0.048 

0.951 

0.023 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation 


Overlay 

Percent 

equation 

Difference 

(cm) 

33.65 

1% 

][|||[||^^ 

4% 

7% 

12.19 

39% 

32.50 

2% 

28.05 

3% 

5% 

6.42 

60% 

31.39 

1% 

26.16 

1% 

1% 

30.25 

4% 

23.91 

4% 

11.77 

7% 

9% 

5% 

6.48 

IBm; 

26.67 

3% 

20.78 

8% 

11.43 

154% 

29.04 

8% 

24.74 

5% 

(microns) 


(mm) 


ESALS 


MPa  (millions) 


Overlay 

EVERPAVE 


Absolute 

Difference 

(cm) 


0.35 


1.13 


1.85 


7.81 


0.50 


0.95 


1.14 


2.42 


17.25 


27.95 


22.78 


31 

70 

1 

12.50 

13.54 

8% 

1.04 

33 

70 

1 

26.00 

26.73 

3% 

0.73 

602  543 

70 

1 

20.00 

20.51 

3% 

0.51 

70 

1 

9% 

0.75 

37 

70 

1 

10% 

2.31 

38 

511  594 

70 

1 

17.00 

17.71 

4% 

39 

291  727 

70 

1 

3.00 

4.19 

40% 

1.19 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Case  Do 
No. 

(microns) 


801 


659 


497 


752 


566 


384 


696 


A 

Esg 

ESALS 

Overlay 

EVERPAVE 

Overlay 

equation 

Percent 

Difference 

Absolute 

Difference 

(mm) 

MPa 

(millions) 

(cm) 

(cm) 

(cm) 

131 


27.00 


25.50 


21.50 


13.00| 


25.50 


20.50 


10.50 


24.00 


18.50 


1.00 


26.00 


24.00 


21.50 


24.50 


20.50 


12.00 


24.50 


23.00 


17.00 


5.50 


20.50 


14.00 


34.501 


33.50 


33.50 


27.56 


26.43 


22.02 


15.07 


23.77 


17.34 


6.81 


22.09 


14.49 


23.91 


35.60 


34.62 


34.81 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Case 

Do 

A 

Esg 

No. 

(microns) 

(mm) 

MPa 

86 

1527 

462 

35 

1591 

450 

35 

1447 

476 

_ 

35 

89 

1529 

35 

90 

1352 

■BB 

35 

91 

1202 

70 

92 

1157 

371 

70 

93 

1174 

367 

70 

94 

377 

70 

95 

370 

70 

96 

386 

70 

97 

1112 

70 

98 

1001 

70 

99 

1064 

70 

100 

925 

HIKQ] 

70 

101 

300 

140 

102 

836 

305 

■H 

103 

847 

303 

■EQI 

ligl 

805 

311 

140 

830 

140 

106 

767 

140 

107 

806 

310 

140 

108 

332 

140 

109 

766 

318 

140 

no 

352 

140 

111 

660 

247 

280 

112 

HES 

250 

280 

113 

IHIQg 

249 

280 

114 

HEEI 

255 

280 

115 

HHE 

251 

280 

116 

585 

261 

280 

HE 

616 

254 

280 

mm 

541 

273 

280 

119 

262 

280 

120 

291 

280 

Overlay 

Percent 

Absolute 

equation 

Difference 

Difference 

(cm) 

(cm) 

32.73 

1% 

0.23 

33.32 

4% 

1.32 

31.81 

4% 

32.62 

9% 

2.62 

30.82 

10% 

2.82 

31.71 

4% 

1.21 

5% 

IHIQ 

7% 

_ m 

8% 

2% 

0.60 

29.31 

3% 

0.81 

29.99 

2% 

0.49 

28.37 

0.87 

29.26 

1.76 

27.27 

9% 

2.27 

28.95 

0% 

0.05 

1% 

0.33 

6% 

1.54 

27.58 

4% 

1.08 

0% 

1% 

27.50 

2% 

0.50 

25.63 

1% 

0.37 

ESALS 


Overlay 

EVERPAVE 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Overlay 

Overlay 

EVERPAVE 

equation 

(cm) 

(cm) 

42.00 

43.57 

39.50 

39.80 

34.00 

33.07 

29.50 

22.11 

41.00 

42.43 

37.50 

37.97 

29.50 

29.79 

18.50 

16.34 

40.00 

41.32 

35.00 

36.08 

25.00 

26.11 

10.19 

32.00 

21.69 

38.86 

30.96 

16.41 

39.00 

40.01 

36.00 

36.59 

32.50 

30.70 

29.50 

21.35 

37.00 

33.00 

25.00 

BQKQI 

19.00 

■QBI 

37.50 

37.87 

32.50 

32.70 

22.50 

23.46 

2.50 

8.82 

36.50 

36.66 

30.00 

30.43 

19.17 

35.23 

26.50 

27.63 

Case  Do 
No. 

(microns) 


1589 


1303 


961 


674 


(mm) 


455 


MPa 


35 


35 


35 


35 


35 


35 


35 


35 


35 


ESALS 


(millions) 


11 

683 

12 

398 

13 

1308 

14 

912 

615 

tm 

553 

17 

1195 

mugg 

18 

787 

19 

453 

ebb 

21 

1105 

22 

907 

423 

23 

675 

486 

24 

494 

517 

25 

1042 

HUgQ 

26 

797 

ebb 

27 

543 

551 

28 

605 

29 

973 

407 

30 

697 

497 

HI 

442 

614 

32 

268 

683 

33 

893 

430 

34 

602 

543 

35 

360 

674 

37 

804 

461 

38 

511 

594 

25% 


3% 


1% 


1% 


12% 


3% 


3% 


4% 


104% 


3% 


11% 


9% 


3% 


3% 


2% 


6% 


28% 


5% 


5% 


9% 


20% 


1% 


1% 


4% 


253% 


0% 


1%| 


4% 


7% 


4% 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Overlay 


Case  Do 
No. 

(microns) 


291 


(mm) 


727 


ESALS 


MPa  (millions) 


Overlay 

EVERPAVE 


Percent 

Absolute 

Difference 

Difference 

(cm) 

1% 

0.12 

1% 

0.48 

8% 

2.41 

26% 

7.44 

4% 

1.35 

3% 

0.95 

11% 

2.50 

163% 

8.97 

1% 

0.38 

3% 

0.79 

3% 

0.54 

1454% 

1% 

3% 

0.74 

1% 

0.23 

3% 

1.02 

0% 

3% 

5% 

3% 

4% 

1.21 

18% 

5.35 

8% 

2.68 

6% 

1.82 

10% 

2.14 

332% 

11.64 

2% 

0.61 

2% 

0.48 

6% 

1.08 

1% 

0.27 

3% 

0.84 

3% 

0.50 

4% 

1.17 

4% 

15.14 


34.11 


28.02 


19.58 


32.27 


25.16 


15.00 


30.17 


22.05 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Case 

Do 

A 

Esg 

ESALS 

Overlay 

Overlay 

Percent 

Absolute 

No. 

EVERPAVE 

equation 

Difference 

Difference 

(microns) 

(mm) 

MPa 

(millions) 

(cm) 

(cm) 

(cm) 

79 

129 

603 

280 

5 

9.50 

10.07 

6% 

0.57 

81 

1749 

435 

35 

5 

43.50 

45.53 

5% 

2.03 

82 

1670 

444 

35 

5 

43.00 

44.55 

4% 

1.55 

83 

1694 

440 

35 

5 

43.00 

4% 

1.73 

84 

1599 

452 

35 

5 

42.00 

4% 

1.58 

85 

1643 

444 

35 

5 

42.50 

43.96 

3% 

1.46 

86 

1527 

462 

35 

5 

41.50 

.  _ _ _ .1 

42.65 

3% 

1.15 

87 

1591 

450 

35 

5 

43.00 

43.25 

1% 

0.25 

40.50 


41.00 


37.50 


41.00 


41.73 


5 

40.00 

5 

39.00 

5 

5 

39.00 

5 

40.50 

5 

38.00 

5 

38.50 

5 

35.50 

5 

38.50 

5 

38.00 

5 

37.00 

•  5 

35.50 

5 

38.00 

5 

36.50 

5 

38.00 

5 

36.00 

5 

36.00 

5 

33.00 

5 

37.00 

5 

36.00 

5 

35.00 

5 

34.00 

35.55 


36.58 


34.22 


38.40 


37.89 


37.99 


36.99 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Case  Do 
No. 

(microns) 


636 


585 


616 


541 


582 


485 


1589 


1303 


961 


674 


1500 


(mm) 


251 


261 


254 


273 


262 


291 


MPa 


280 


280 


280 


280 


280 


280 


ESALS  Overlay 

EVERPAVE 

(millions)  (cm) 


36.00 


35.00 


36.00 


34.00 


34.00 


31.00 


Overlay 

equation 

(cm) 


37.55 


35.93 


36.87 


34.54 


35.78 


32.81 


Percent 

Difference 


Absolute 

Difference 

(cm) 


1.55 


0.93 


0.87 


0.54 


1.78 


1.81 


23 

675 

486 

70 

10 

7% 

2.52 

24 

494 

517 

70 

10 

25% 

8.37 

25 

1042 

391 

70 

10 

42.00 

43.24 

3% 

1.24 

797 

10 

38.00 

38.94 

2% 

0.94 

543 

10 

30.00 

31.45 

1.45 

28 

358 

605 

70 

10 

23.50 

19.40 

17% 

4.10 

29 

973 

407 

70 

10 

2% 

0.85 

30 

697 

497 

70 

10 

3% 

1.02 

31 

442 

614 

70 

10 

1% 

0.26 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


(microns) 


268 


893 


602 


360 


804 


511 


291 


801 


659 


A 

Esc 

ESALS 

Overlay 

EVERPAVE 

Overlay 

equation 

Percent 

Difference 

Absolute 

Difference 

(mm) 

MPa 

(millions) 

(cm) 

(cm) 

(cm) 

43 

403 

140 

10 

36.00 

7% 

2.63 

44 

429 

140 

10 

33.50 

23% 

7.66 

45 

752 

140 

10 

40.63 

3% 

1.13 

46 

566 

140 

10 

36.22 

2% 

0.72 

■B 

384 

470 

140 

10 

^/ggUfj^ 

29.27 

6% 

1.77 

48 

262 

140 

10 

22% 

5.25 

49 

696 

140 

10 

Hugos 

3% 

1.11 

50 

140 

10 

33.98 

4% 

1.52 

286 

392 

10 

26.50 

28.41 

7% 

1.91 

68 

203 

435 

280 

10 

24.00 

gggg^ 

4.59 

69 

526 

277 

280 

10 

38.00 

gggg^ 

0.38 

70 

354 

352 

280 

10 

33.00 

2% 

0.70 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


(microns) 


(mm) 


Esg  ESALS 
MPa  (millions) 


Overlay 

EVERPAVE 

(cm) 


Overlay  Percent 
equation  Difference 

(cm) 


45.00 


44.00 


42.50 


45.50 


44.00 


46.00 


43.50 


44.50 


40.50 


43.50 


42.50 


42.00 


40.00 


43.00 


1.50 

3.501 


41.00 


45.13 


45.36 


44.34 


Absolute 

Difference 

(cm) 


0.85 


0.45 


1.06 


0.23 


0.95 


87 

1591 

450 

35 

10 

49.00 

3% 

1.48 

88 

^bee& 

476 

35 

46.50 

1% 

0.49 

89 

460 

35 

HBES 

47.50 

46.82 

1% 

0.68 

90 

497 

35 

10 

43.50 

3% 

1.52 

91 

_ 1 

364 

70 

10 

46.00 

0% 

0.09 
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Table  A-3  Comparison  of  Overlay  Required:  EVERPAVE  vs  Equation,  Cont. 


Case 

No. 

Do 

(microns) 

A 

(mm) 

Esg 

MPa 

ESALS 

(millions) 

Overlay 

EVERPAVE 

(cm) 

Overlay 

equation 

(cm) 

Percent 

Difference 

Absolute 

Difference 

(cm) 

109 

766 

318 

140 

10 

41.50 

40.86 

2% 

0.64 

653 

10 

38.00 

38.50 

1% 

0.50 

■m 

660 

mg 

10 

42.68 

3% 

1.18 

1F3 

643 

250 

280 

10 

40.50 

42.16 

4% 

1.66 

113 

649 

249 

280 

10 

39.50 

7% 

114 

617 

255 

280 

10 

38.00 

imD^ 

9% 

115 

636 

280 

10 

41.83 

3% 

1.33 

116 

585 

1^^ 

280 

10 

40.20 

2% 

0.70 

117 

616 

254 

280 

10 

41.14 

0% 

0.14 

118 

541 

273 

280 

10 

38.50 

38.81 

0.31 

119 

262 

280 

10 

39.00 

40.06 

3% 

1.06 

120 

485 

291 

280 

10 

36.00 

37.09 

3% 

1.09 

Total 

4310% 

461.71 

Average 

13% 

1.411971 

APPENDIX  B 
GRAPHS 


B.1  INTRODUCTION 

The  best  way  to  see  if  there  is  a  relationship  between  two  parameters  is  to  graph 
one  versus  the  other  on  an  x-y  scatter  plot.  A  regression  equation  defines  the  “best-fit” 
straight  line  that  represents  the  relationship  between  the  two  variables.  Often,  variables 
relate  to  each  other  in  a  non-linear  fashion.  In  this  case,  transformations  can  be  made  to 
one  or  both  of  the  variables  in  order  to  create  a  more  linear  relationship  Figure  B-1 
shows  four  typical  relationships  between  two  parameters.  For  each  of  these  examples,  a 
generic  equation  can  be  written.  These  relationships  and  their  equations  are: 

a.  Linear  y  =  bo  +  biX 

b.  Exponential  y  =  bobi* 

c.  Power  y  =  boX*’' 

d.  Hyperbolic  y  =  bo  +  b 

where  bo  and  bi  are  adjusted  to  best  represent  the  relationship  between  x  and  y. 

B.2  GRAPHS 

Figures  B-2  through  B-1 1  show  the  graphs  produced  during  the  first  stage  of 
research  for  this  thesis.  These  graphs  relate  different  measured  pavement  responses  such 
as  Do  and  Area  Parameter  to  the  critical  strains,  Sj  8v. 
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Linear  Relationship  Between  x 


and  y 


Power  Relationship  Between  x 


and  y 


Exponential  Relationship 


Between  x  and  y 


Hyperbolic  Relationship 
Between  x  and  y 


Figure  B-1  Typical  Relationships  Between  Variables  (after  Mahoney,  1994) 
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Subgrade  Modulus  =  35  MPa 


Area  Parameter  (mm) 


Figure  B-3  Area  Parameter  vs  for  Subgrades  with  Modulus  =  35  MPa 


Subgrade  Modulus  =  70  MPa 


Area  Parameter  (mm) 


Figure  B-4  Area  Parameter  vs  8v  for  Subgrades  with  Modulus  =  70  MPa 
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Subgrade  Modulus  =  140  MPa 
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Figure  B-5  Area  Parameter  vs  8y  for  Subgrades  with  Modulus  =  140  MPa 


Subgrade  Modulus  =  280  MPa 


Area  Parameter  (mm) 


Figure  B-6  Area  Parameter  vs  8v  for  Subgrades  with  Modulus  =  280  MPa 
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Subgrade  Modulus  =  35  MPa 
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Figure  B-8  Area  Parameter  vs  for  Subgrades  with  Modulus  =  35  MPa 
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Figure  B-9  Area  Parameter  vs  8*  for  Subgrades  with  Modulus  =  70  MPa 
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Subgrade  Modulus  =  140  MPa 
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Figure  B-10  Area  Parameter  vs  St  for  Subgrades  with  Modulus  =  140  MPa 
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Figure  B-11  Area  Parameter  vs  e*  for  Subgrades  with  Modulus  =  280  MPa 


APPENDIX  C 


REGRESSION  ANALYSIS 

C.l  INTRODUCTION 

Regression  analysis  is  a  powerful  tool  used  to  establish  relationships  that  exist  in 
a  set  of  data.  Regression  analysis  is  used  extensively  in  this  thesis  therefore,  this 
Appendix  is  included  to  provide  a  cohesive  explanation  of  the  regression  analysis 
proceedure.  Section  3.2.2  provides  a  brief  description  of  regression  analysis.  This 
section  will  go  into  greater  detail. 

C.2  CORRELATION 

Correlation  is  a  method  of  quantifying  the  association  of  two  variables.  The 
Pearson  product  moment  correlation  coefficient,  r,  is  commonly  used.  The  basic 
equation  for  r  is  (Ryan  et  al.,  1985): 

_  ^(x-r)(y-y) 

while  the  formula  normally  used  to  compute  r  is  (Ryan  et  al.,  1985); 

Sv — — 

Its  value  is  always  between  -1  and  +1.  On  an  x-y  plot  of  the  data,  a  positive  correlation 
indicates  y  tends  to  increase  as  x  increases,  commonly  known  as  a  positive  slope  in 
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mathematics.  The  converse  is  true  for  a  negative  correlation.  It  indicates  a  negative 
slope:  y  tends  to  decrease  as  x  increases.  If  no  trend  can  be  seen  at  all,  r  will  equal  0; 
although,  an  r  equal  to  0  does  not  always  indicate  a  lack  of  association  between  x  and  y. 

It  indicates  a  lack  of  a  linear  association  between  x  and  y.  A  curved  association  is  an 
example  of  a  relationship  that  cannot  be  identified  by  the  correlation  coefficient. 

Because  of  this  limitation,  a  plot  of  the  data  on  an  x-y  scatter  plot  provides  an  excellent 
picture  of  any  association  that  exists  between  the  variables.  The  correlation  coefficient 
should  not  be  used  exclusively  to  determine  if  two  variables  are  related  by  some 
function. 

C.3  REGRESSION 

Whereas  correlation  provides  a  quantitative  measure  of  relationship  between  two 
variables,  regression  provides  the  equation  to  describe  the  relationship  mathematically. 
Simple  regression  involves  only  two  variables  and  is  the  most  straightforward  to  explain. 
Simple  regression  calculates  the  “best-fit”  line  on  an  x-y  plot  of  the  two  variables. 

Figure  C-1  illustrates  the  “best-fit”  line  using  an  example  of  height  vs  weight  for  5 
people. 

The  method  of  least  squares  is  used  to  calculate  the  “best-fit”  line.  The  equation 
for  any  straight  line  follows  the  form  y  =  a  +  hx,  where  x  and  are  the  variables  and  a 
and  b  are  the  constants  defining  the  relationship.  Using  the  least  squares  method,  a  and  h 
are  calculated  by  (Ryan  et  al.,  1985): 
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Height  vs  Weight 


Figure  C-1  Illustration  of  the  “Best  Fit”  Line 

Multiple  regression  involves  more  than  two  variables.  The  development  of  a 
regression  equation  is  more  rigorous  using  multiple  variables;  therefore,  statistical 
computer  programs  such  as  Minitab  are  used  to  calculate  the  equation.  Most  of  the 
equations  developed  in  this  thesis  utilized  multiple  regression,  they  can  be  found  in 
Appendix  D. 
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The  parameters  used  to  quantify  the  usefulness  of  a  regression  equation  are 
explained  in  Chapter  3,  Section  3.2.2.  These  are  R^,  RMSE  and  N. 


APPENDIX  D 


MINITAB  OUTPUT 


Table  D-1  Regression  Equations  for  et  and  log 


Equation 

8t  =  984-  1.13  A  -  l.lOEsG 
8t  =  -79.9  +  0.434  Do  +  0.373  Ego 
log  St  =  -4.23  +  1.80  log  Do  +  0.898  log  Ego 
log  St  =  -2.3 1  +  1.60  log  Do  +  0.00273  Eg^ 
log  St  =  5.91  -  0.00353  A  -  0.877  log  Ego 
log  St  =  4.59 -0.00351  A -0.00335  Ego 
log  St  =  18.0 -5.12  log  A -0.863  log  Ego 
log  St  =  16.6  -  5.10  log  A  -  0.00330  Egg 
log  St  =  3.88 -0.00290  A 


Wa 

97^ 

90 

80 

80 

STA 

06 

OJ 


RMSE 

80.51 

48.46 

0.06465 

0.1080 

0.1358 

0.1593 

0.1470 

0.1678 

0.3081 


N 

”41 

O 

IT 

"TT 

O 

IT 

”41 

IT 


Table  D-2  Regression  Equations  for  log  Sy 


Equation 

mn 

RMSE 

N 

log  Sv  =  7.04  -  0.00339A  -  1.34  log  Ego 

94.4 

0.01148 

120 

log  Sv  =  15.3  -  3.67  log  A  -  1.42  log  Ego 

91.0 

0.1450 

120 

log  Sv  =  12.7  -  3.52  log  A  -  0.00485  Ego 

86.7 

log  Sv  =  4.88  -  0.00323A  -  0.00451EgG 

88.8 

0.1616 

120 

log  Sv  =  1.89  +  0.00112  Do  -  0.000833  EgG 

81.6 

0.2075 

120 

log  Sv  =  6.0  -  0.00284  A  -  1.03  log  EgG  +0.000233  Dq 

120 

0.1100 

120 

log  Sv  =  8.73  -  0.150  -  1.38  log  EgG 

120 

0.1247 

120 

Table  D-3  Regression  Equations  for  8, 


Equation 

RMSE 

n 

s 

^  =  3964-4.3  A -7.54  Ego 

69.5 

450.2 

8 

^  =  5994  -  1.88  (Af"  -  7.74  Ego 

68.2 

460.0 

120 

8 

v  =  -511  +  1.97  Do +  0.788  Ego 

96.0 

163.9 

120 

8 

^  =  -88.5  +  1.83  Do -0.472  A 

96.1 

161.6 

120 

8 

V  =  -204  +  1.87  Do  -  0.359  A  +  0.252  Ego 

96.1 

161.9 

120 

8 

^--337  +  1.88  Do 

95.4 

174.5 

120 

Table  D-4  Regression  Equations  for  Area  Parameter  Load  Adjustment 


Case 

Number 


50 


Equation 

N 

■ 

RMSE 

A  =  380  +  0.968  (Load) 

7 

98.8 

3.337 

A  =  310 +  0.982  (Load) 

7 

99.2 

A  =  567  +  0.835  (Load) 

7 

95.8 

A  =  217+  1.16  (Load) 

7 

99.9 

1.406 

A  =  201  +  1.23  (Load) 

m 

99.8 

1.527 

A  =  395  +  1.36  (Load) 

m 

97.9 

6.297 

A  =  761  +0.715  (Load) 

7 

93.1 

6.155 

A  =  410+  1.18  (Load) 

7 

98.9 

3.915 

A  =  485  +  1.37  (Load) 

7 

90.3 

14.18 

A  =  468  +  0.888  (Load) 

7 

99.1 

2.674 

Range  of  slope  of  line:  0.715  -  1.37 
Average  slope:  1.069 

Use  average  slope  as  Area  Parameter  adjustment  factor  for  load 
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Table  D-5  Regression  Equations  for  Dq  Load  Adjustment 


Case 

Number 

Equation 

N 

R" 

RMSE 

50 

Do  =  76.8  +  9.93  (Load) 

7 

99.8 

12.95 

70 

Do  =  73.3  +  6.78  (Load) 

7 

99.7 

11.97 

55 

Do  =  20.8  +  5.42  (Load) 

7 

100 

2.672 

65 

Do  =  225  +  8.13  (Load) 

7 

98.7 

30.00 

115 

Do  =  303  +  7.65  (Load) 

7 

98.2 

32.38 

1 

Do  =  316 +  31.0  (Load) 

7 

99.8 

40.87 

40 

Do  =  8.01  +3.76  (Load) 

7 

100 

0.7115 

89 

Do  =  259  +  30.9  (Load) 

7 

99.8 

42.97 

10 

Do  =  85.5  +  23.5  (Load) 

7 

100 

15.35 

13 

Do  =  151  +28.4  (Load) 

7 

99.9 

28.74 

Range  of  slope  of  line;  varies  with  initial  value  of  Dq 
Develop  equation  to  predict  the  slope  of  the  line  based  on  load. 


Table  D-6  Regression  Equations  for  Slope  of  Do  Adjustment  for  Load 


Equation 

N 

RMSE 

slope  =  -0.91  +  0.0369  (Do)2o 

10 

90.2 

3.788 

slope  =  -1.01  +  0.0265  (Do)3o 

10 

93.8 

3.022 

slope  =  -0.98  +  0.0209  (Do)4o 

10 

95.9 

2.463 

slope  =  -0.92  +  0.0174  (Do)5o 

10 

97.2 

2.030 

slope  =  -0.83  +  0.0149  (Do)6o 

10 

97.8 

1.686 

slope  =  -0.737  +  0.0131  (Do)7o 

10 

98.7 

1.402 

slope  =  -0.643  +  0.0117  (Do)8o 

10 

99.1 

1.164 

1/pslope  =  8.51  +  0.970  (Load) 

7 

99.9 

0.5917 

intercept  =  -1.13  +  0.00535  (Load) 

(all  loads) 

7 

75.4 

0.0722 

intercept  =  -1.27  +  0.00758  (Load) 

(loads  30-80  kN) 

7 

97.4 

0.02611 

Use  the  equation  slope  =  intercept  +  psiope  (Load)  to  calculate  the  slope  of  the  Load 
Adjustment  Factor  for  Do 

Where  intercept  =  -1 .27  +  0.00758  (Load)  for  loads  30-80  kN 
intercept  =  -0.91  for  20  kN  load 

and  psiope  =  1/  (8.51  +  0.970  (Load))  for  all  loads 
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Figure  D-7  Regression  Equations  for  Temperature  Shift  Factor  for  A 


Equation 


Asf  =  4.96  + 0.0375  Esg 


Agp  — 11.2  +  10.6  log  Ego 


AgF=  14.5-345(l/EgG) 


AgF  =  7.09  +  0.000107(EgG) 


AgF  =  0.45  +  0.880(EgG) 


AgF  =  -2.77  +  0.880(EgG)  +  0-238  Tag 


AgF  =  6.24  +  0.880(EgGr"  -  33.3(1/Tac) 


AgF  =  -6.64  +  0.880(EgG)  +  2.16(Tac) 


log  AgF  =  0.849  +  0.0361(EgG)  -  1-93(1/Tac) 
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Table  D-8  Regression  Equations  for  Area  Temperature  Shift  Factor  Calculation 


Esg 

Tac 

Regression  Equation 

R^ 

RMSE 

N 

35 

2.5 

1/A  =  0.00219  +  0.000002(Temp) 

97.9 

8.3E''’ 

5 

70 

2.5 

1/A  =  0.00262  +  0.000002(Temp) 

94.3 

5 

140 

2.5 

1/A  =  0.00318  +  0.000003(Temp) 

92.2 

2.4F^ 

5 

280 

2.5 

1/A  =  0.00387  +  0.000003(Temp) 

90.5 

3.1E'^ 

5 

35 

3.75 

1/A  =  0.00205  +  0.000004(Temp) 

97.2 

1.8E‘^ 

5 

70 

3.75 

1/A  =  0.00242  +  0.000005(Temp) 

96.4 

2.8E’^ 

5 

140 

3.75 

1/A  =  0.00291  +  0.000007(Temp) 

95.7 

4.2E'^ 

5 

280 

3.75 

1/A  =  0.00353  +  0.000009(Temp) 

94.9 

5.8E'^ 

5 

35 

5.0 

1/A  =  0.00191  +  0.000005(Temp) 

98.6 

1.8E'^ 

5 

70 

5.0 

1/A  =  0.00223  +  0.000007(Temp) 

98.3 

2.7F^ 

5 

140 

5.0 

1/A  =  0.00265  +  0.0000  lO(Temp) 

98.1 

4.0F^ 

5 

280 

5.0 

1/A  =  0.00319  +  0.000013(Temp) 

97.8 

5.7F^ 

5 

35 

10.0 

1/A  =  0.00154  +  0.000008(Temp) 

100.0 

3JF^ 

5 

70 

10.0 

1/A  =  0.00172  +  0.00001  l(Temp) 

99.9 

7.0F'’ 

5 

140 

10.0 

1/A  =  0.00196  +  0.0000 15(Temp) 

99.9 

I.IF^ 

5 

280 

10.0 

1/A  =  0.00229  +  0.000021(Temp) 

100.0 

1.3F^ 

5 

35 

20.0 

1/A  =  0.00129  +  0.000007(Temp) 

97.8 

3.2F^ 

5 

70 

20.0 

1/A  =  0.00138  +  0.00001  l(Temp) 

97.6 

4.8F^ 

5 

140 

20.0 

1/A  =  0.00150  +  0.000016(Temp) 

97.6 

7.0F' 

5 

280 

20.0 

1/A  =  0.00167  +  0.000022(Temp) 

97.9 

9.2F^ 

5 

35 

40.0 

1/A  =  0.001 19  +  0.000007(Temp) 

94.6 

4.9F^ 

5 

70 

40.0 

1/A  =  0.00124  +  0.00001  l(Temp) 

94.7 

7.2F^ 

5 

140 

40.0 

1/A  =  0.00132  +  0.000016(Temp) 

95.1 

1.0F‘' 

5 

280 

40.0 

1/A  =  0.00144  +  0.000022(Temp) 

96.0 

1.3F‘* 

5 
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Table  D-9  Regression  Equations  for  D©  Temperature  Shift  Factor  Calculation 


Table  D-10  Regression  Equations  for  Overlay,  N  =  327 
(without  zero  values  for  overlay  required) 
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Equation 

N 

RMSE 

overlay  =  -135  +  451ogDo  +  14.9logEsG  +  14.21ogESALs 

327 

93.7 

2.686 

overlay  =  -541  +  84.91ogDo  +  82.41ogA  +  54.31ogEsG  + 
14.51ogESALs 

327 

95.7 

2.239 

overlay  =  111  -  0.084A  -  26.31ogEsG  13.81ogESALs 

327 

90.7 

3.273 

overlay  =  -167  +  50.81ogDo  +  0.01 13A  +  20.41ogEgG  + 
14.21ogESALs 

327 

93.8 

2.682 

Table  D-11  Regression  Equations  for  Overlay,  N  =  360 
(including  zero  values  for  overlay  required) 


Equation 

N 

H 

RMSE 

overlay  =  -8.06  +  0.0291Do  +  0.0501Esg  +  1.53ESALs 

360 

77.2 

6.42 

overlay  =  -137  +  45.41ogDo  +  15.41ogEsG  +  13.81ogESALs 

360 

93.8 

3.337 

log(overlay)  =  -2.53  +  1.061ogDo  +  0.4131ogEgG  + 
0.2311ogESALs 

360 

66.5 

0.1591 

overlay  =  -109  +  44.  llogDo  +  0.0505Esg  +  1.53ESALs 

360 

91.6 

3.906 

overlay  =  -429  +  72.61ogDo  +  61.41ogA  +  43.21ogEgG  + 
13.81ogESALs 

360 

94.9 

3.032 

overlay  =  69.8  -  0.0846A  -  0.0893Esg  +  1.53ESALs 

360 

84.9 

5.224 

overlay  =  115  -  0.089A  -  27.41ogEsG  +  1.53ESALs 

360 

89.9 

4.266 

log(overlay)  =  3.24  -  0.00 194 A  -  0.5511ogEsG  + 
0.0241ogESALs 

360 

61.0 

0.1717 

overlay  =  83  -  0.0874A  -  2.2(Esg)”  ^  +  1.53ESALs 

360 

88.1 

4.641 

overlay  =  61.6  -  0.0788A  -  0.0002(Esg)^  +  1.53ESALs 

360 

78.8 

6.192 

overlay  =  -6.82  +  16410(1/A)  -  0.102Esg  +  1.53ESALs 

360 

75.9 

6.602 

overlay  =  274  -  91.61ogA  -  0.0976Esg  +  1.53ESALs 

360 

82.2 

5.684 

overlay  =  48.4  -  0.00008(A)^  -  0.0798Esg  +  1.53ESALs 

360 

84.5 

5.302 
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Table  D-11  Regression  Equations  for  Overlay,  N  =  360 
(including  zero  values  for  overlay  required),  cont.. 


Equation 

N 

lai 

RMSE 

overlay  =110-  3.73A”‘^  -  0.0937EgG  +  1.53ESALs 

360 

84.0 

5.384 

overlay  =  46.6  -  0.0863 A  +  918(1/Esg)  +  1.53ESALs 

360 

87.6 

4.730 

overlay  =  159  -  3.91A“^  -  28.51ogEgG  +  1.53ESALs 

360 

88.8 

4.506 

overlay  =  331  -  95.71ogA  -  29.41ogEsG  +  1.53ESALs 

360 

86.5 

4.944 

overlay  =  38.2  +  16834(1/A)  -  29.81ogEsG  +  1.53ESALs 

360 

78.8 

6.200 

overlay  =  88.3  -  0.00008A"-^  -  24.81ogEsG  +  1.53ESALs 

360 

89.4 

4.377 

(overlay)^  =  4810  -  3.83A  -  12951ogEsG  +  87.1ESALs 

360 

89.2 

207.4 

overlay  =  115-  0.0890A  -  27.41ogEsG  +  13.81ogESALs 

360 

91.0 

4.045 

overlay  =  129  -  0.0890A  -  27.41ogEsG  -  14.0(l/ESALs) 

360 

89.9 

4.274 

overlay  =  109  -  0.0890A  -27.41ogEsG  +  6.52(ESALs)'^  ^ 

360 

90.8 

4.071 

overlay  =  118-  0.0890A  -27.41ogEsG  +  0. 123(ESALs)'^ 

360 

87.2 

4.808 

APPENDIX  E 


NONDESTRUCTIVE  TESTING  DEVICES 

E.l  INTRODUCTION 

Many  summaries  of  payment  nondestructive  testing  devices  have  been  written. 
This  appendix  is  included  as  a  general  overview  of  the  most  commonly  used  devices; 
past,  present,  and  future.  As  a  reminder,  only  devices  that  can  measure  the  deflection 
basin  can  be  used  to  calculate  the  Area  Parameter  used  frequently  in  this  thesis.  Most 
devices  are  able  to  measure  the  maximum  deflection  at  the  center  of  the  load  (Dq). 
Chapter  2  includes  a  summary  of  which  devices  are  most  commonly  used  in  the  United 
States  (Table  2-1). 

E.2  STATIC  DEFLECTION  EQUIPMENT 

Devices  that  measure  the  pavement’s  response  to  a  slowly  applied  load  are 
classified  as  static  deflection  equipment.  The  Benkelman  Beam  was  the  most  widely 
used  static  deflection  device  resulting  in  much  of  the  earlier  deflection  based  overlay 
designs  being  based  on  Benkelman  Beam  data  (Smith  and  Lytton,  1984). 

E.2.1  Benkelman  Beam 

The  Benkelman  Beam  consists  of  a  3.66  meter  (12  ft)  beam  with  a  probe  on  one 
end  that  rests  on  the  pavement  between  the  rear  dual  tires  of  a  loaded  truck  axle.  The 
other  end  is  supported  by  stationary  legs  that  are  ideally  outside  the  influence  of  the  truck 
load.  There  is  a  pivot  point  2.44  meters  (8  ft)  from  the  probe  and  the  remaining  1.22 
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meters  (4  ft)  of  the  beam  activates  a  dial  gauge  for  measuring  the  displacement.  The 
displacement  is  read  from  the  dial  as  the  truck  is  slowly  driven  away  from  the  probe. 

The  maximum  dial  reading  is  recorded  as  the  rebound  deflection  (Newcomb,  1986,  p.  8). 

E.2.2  Plate  Bearing  Test 

The  plate  bearing  test  loads  a  circular  plate  on  the  pavement  using  a  hydraulic 
jacking  system  that  reacts  against  the  frame  of  a  truck.  The  plate’s  deflection  is 
measured  by  a  dial  gauge  which  is  mounted  on  a  stand  placed  as  far  away  from  the  plate 
as  possible.  Curvature  in  the  plates  is  reduced  by  using  a  stack  of  successively  smaller 
diameter  steel  plates  (Smith  and  Lytton,  1984). 

E.3  AUTOMATED  BEAM  DEFLECTION  EQUIPMENT 

This  category  of  deflection  equipment  automates  the  Benkelman  Beam  process. 
E.3.1  La  Croix  Deflectograoh 

The  La  Croix  Deflectograph  consists  of  a  truck  with  a  pivot/beam  assembly 
mounted  on  it.  The  length  and  geometry  of  the  beam  vary  with  different  wheel  bases  for 
different  trucks.  The  measurement  process  begins  with  the  beam  being  placed  in  front  of 
the  wheel.  As  the  truck’s  wheel  approaches  the  beam  tip,  the  beam  rotates  about  the 
pivot.  The  rotation  is  measured  until  the  wheel  has  gone  approximately  21.9  cm  (0.72  ft) 
past  the  beam  tip.  The  beams  and  frame  are  then  repositioned  to  be  in  front  of  the  wheel 
again  to  begin  another  measurement.  Measurements  can  be  made  every  3.5  to  6  meters 
(1 1.4  to  19.7  ft)  depending  on  the  particular  truck  (Smith  and  Lytton,  1984).  This  device 
was  used  extensively  in  Europe  but  did  not  gain  popularity  in  the  United  States. 
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E.3.2  California  Traveling  Deflectometer 

The  California  Traveling  Deflectometer  is  very  similar  to  the  La  Croix  Deflectograph 
and  has  been  used  extensively  by  the  California  Department  of  Transportation.  It 
basically  consists  of  two  Benkelman  Beam  probes  mounted  on  a  semi-trailer.  The 
measurement  is  conducted  in  much  the  same  manner  as  the  La  Croix  Deflectograph.. 

E.4  STEADY  STATE  DYNAMIC  DEFLECTION  EQUIPMENT 

This  class  of  deflection  equipment  uses  dynamic  force  generators  to  produce  a 
sinusoidal  vibration  in  the  pavement.  The  most  commonly  used  models  are  the 
Dynaflect  and  the  Road  Rater. 

E.4.1  Dvnaflect 

The  Dynaflect  transmits  a  load  to  the  pavement  through  two  16  inch  diameter 
rubber-coated  steel  wheels,  10  centimeters  (4  in)  wide,  spaced  50.8  centimeters  (20  in) 
center  to  center.  The  dynamic  force  is  created  by  counter-rotating,  eccentric  flywheels 
which  generate  a  4448  newton  (1,000  pound)  peak-to-peak  dynamic  force  which  is  offset 
by  the  machine’s  static  weight  of  approximately  746.5  kilograms  (2,000  pounds). 

Surface  deflections  are  measured  by  five  geophones  spaced  at  30.5  centimeter  (1  ft) 
intervals  from  the  center  of  the  load.  Deflection  measurements  are  recorded  by  a 
computer  located  inside  the  tow  vehicle.  The  vibration  frequency  and  load  cannot  be 
changed  (Newcomb,  1986)  (Smith  and  L5qton,  1984). 
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E.4.2  Road  Rater 

The  Road  Rater,  like  the  Dynaflect,  is  a  steady  state  dynamic  device.  The  load  is 
applied  to  the  pavement  surface  through  a  steel  loading  plate  which  is  approximately  10 
centimeters  by  18  centimeters  (4  in  by  7  in).  Three  models  are  available  from 
Foundation  Mechanics,  Inc.,  the  400B  with  a  maximum  static  load  of  10,676  newtons 
(2,400  pounds),  the  model  2000  with  a  maximum  static  load  of  16,903  newtons  (3,800 
pounds)  and  the  model  2008  with  a  maximum  static  load  of  25,800  newtons  (5,800 
pounds).  The  static  load  is  created  by  a  hydraulic  system  acting  on  the  trailer  weight 
while  the  dynamic  force  is  produced  by  a  hydraulic  actuator  oscillating  a  steel,  lead-filled 
mass.  Each  Road  Rater  model  has  the  capability  of  adjusting  the  loads  and  the 
frequency.  Four  geophones  spaced  at  30.5  centimeter  (1  ft)  intervals  from  the  center  of 
the  loading  plate  measure  the  pavement  deflections  (Smith  and  Lytton,  1984)  (Newcomb, 
1986). 

E.5  IMPULSE  DEFLECTION  EQUIPMENT 

This  category  of  equipment  contains  all  Falling  Weight  Deflectometers,  different 
models  of  which  are  manufactured  by  Dynatest,  KUAB,  and  Phoenix.  These  machines 
all  deliver  a  transient  force  impulse  to  the  pavement,  simulating  the  passage  of  a  truck 
wheel  over  the  pavement  section. 

E.5.1  Falling  Weight  Deflectometer 

The  FWD  is  a  very  precise  measurement  system  which  has  been  shown  to  closely 
model  the  effects  of  truck  loading  on  pavements.  A  heavy  load  is  dropped  from  a 
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predetermined  height  onto  a  rubber  buffer  system  creating  a  load  pulse  in  approximately 
a  half-sine  wave  form  of  25  to  30  millisecond  duration.  The  load  in  transmitted  to  the 
pavement  through  a  load  distributing  steel  plate  situated  directly  on  the  pavement.  The 
plate  typically  measured  30  centimeters  (11.8  in)  in  diameter,  however  other  diameter 
plates  are  also  available  (Smith  and  L54ton,  1984).  An  array  of  geophones  or 
accelerometers  is  placed  from  the  center  of  the  load  plate  radially  outward  along  the  load 
axis.  Depending  on  whether  the  sensors  measure  the  velocity  or  acceleration  of  the 
impulse  waves,  the  data  is  integrated  once  or  twice  respectively  to  compute  the 
deflections  due  to  the  loading.  The  array  of  sensors  capture  data  which  provides  a 
picture  of  the  deflection  basin  due  to  the  load.  Figure  E-1  shows  the  layout  of  a  typical 


FWD  configuration. 
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Loading  Wheel  Contact  Area 


Geophone  (Deflection  Sensor) 


Figure  E-1:  Typical  FWD  Configuration  of  Loading  Plate  and  Geophones 

(WSDOT  Guide,  1995) 

Figure  E-2  illustrates  a  deflection  basin  measured  by  a  FWD.  This  figure  also  shows 
how  different  layers  affect  different  parts  of  the  deflection  basin,  i.e.  the  surface  course 
stiffness  is  most  related  to  Dq  while  the  subbase  quality  is  reflected  in  D3. 
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Figure  E-2  Dlustration  of  FWD  DeHection  Basin  (after  WSDOT  Guide,  1995) 

If  the  structure  of  the  pavement  is  known  (thickness  and  material  in  each  layer)  the 
elastic  modulus  of  each  layer  can  be  estimated  through  the  iterative  process  of 
backcalculation  which  will  be  discussed  in  greater  detail  later  in  this  chapter. 
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E.6  LASER  ROLLING  WEIGHT  DEFLECTOMETER 

A  rolling  wheel  pavement  deflectometer  was  first  reported  and  tested  in  the  mid 
1970’s  (Johnson,  1995).  Harr  began  work  in  the  1970’s  with  several  graduate  students 
evaluating  the  use  of  the  deflection  basin  to  evaluate  airfield  pavements.  They  began  the 
process  of  developing  a  system  that  could  measure  pavement  deflections  from  a 
noncontact  sensor  mounted  on  a  beam  moving  above  the  pavement.  A  rigid  beam  was 
constructed  with  light-emitting  diodes  (LEDs),  however,  ambient  light  was  found  to 
adversely  affect  the  performance  of  the  LEDs  so  they  were  subsequently  replaced  with 
lasers  (Bush,  et  al.,  1983).  In  1984  Bush,  Cox  and  Hall  researched  a  system  of  four 
sensors  on  a  rigid  beam.  Other  efforts  by  Harr  (1975)  and  Bush  et  al.  (1983)  used  the 
four  sensors  on  a  beam  approach,  however,  critical  errors  due  to  beam  deflections  were 
unresolved  (Johnson,  1995).  What  was  needed  was  a  beam  long  enough  to  establish  a 
reference  elevation  of  the  pavement  outside  the  deflection  basin,  yet  rigid  enough  not  to 
create  error  due  to  the  deflection  in  the  beam  itself  (or  some  method  of  calculating  the 
beam  deflection  at  any  given  time  and  adjusting  the  pavement  deflection  accordingly). 
E.d.l  Two-Point  Deflection  Measurement 

The  principle  for  measuring  pavement  deflections  using  a  beam  passing  over  the 
pavement  is  to  determine  two  heights  from  a  floating  reference  datum;  one  on  the 
unloaded  pavement  and  the  second  at  the  same  point  on  the  pavement  at  the  time  the 
loaded  wheel  passes  (Johnson,  1995).  The  deflection  is  then  the  difference  between  these 
two  calculated  heights.  The  Transportation  and  Road  Research  Laboratory  (TRRL)  in 
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Crowthome,  England,  conducted  research  on  a  high-speed  profilometer  using  4  sensors 
on  a  rigid  beam.  The  solution  to  the  noncontact  measurement  system  used  by  many 
rolling  deflectometers  was  derived  from  this  TRRL  research  (Bush  et  al.,  1983).. 

In  order  to  measure  the  distance  from  the  reference  datum  to  the  undeflected 
pavement,  three  sensors  are  required  outside  the  influence  of  the  applied  load  (beyond 
the  deflection  basin).  Bush,  Hall  and  Harr  (1983)  described  the  profile  algorithm  used  to 
measure  the  undeflected  pavement  profile.  The  description  of  the  algorithm  that  follows 
uses  nomenclature  consistent  with  that  used  by  Johnson  (1995)  in  the  RWD  Final  Report. 
The  undeflected  height  above  the  pavement  is  calculated  from  an  initial  set  of 
measurements  recorded  from  the  three  equally  spaced  lead  sensors,  termed  sensors  A,  B, 
C  as  shown  in  Figure  E-3.  The  height  measured  by  sensors  A,  B  and  C  are  termed  A,  B, 
and  C  respectively.  Using  similar  triangles  and  a  level  reference  datum,  a  virtual  height, 
h,  can  be  calculated  at  sensor  C  by  the  equation: 

A  -  (A-B)  =  C-h  +  (A-B) 
solving  for  h: 


\i=A-2B+C 
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Figure  E-3  First  Point  Measurement  (after  Johnson,  1995) 

A  fourth  sensor,  D  is  located  adjacent  to  the  load  wheel  and  measures  the  deflection  due 
to  the  applied  load.  The  RWD  is  designed  so  that  all  four  sensors  are  set  at  9  ft  spacing. 
The  virtual  deflected  height,  h'  is  calculated  when  sensor  D  is  directly  over  the  spot 
sensor  C  was  at  when  h  was  calculated  as  shown  in  Figure  E-4. 
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Figure  E-4  Second  Point  Measurement  (after  Johnson,  1995) 

Again,  using  similar  triangles  and  a  level  reference  datum  the  equation  to  calculate  h'  is 

B -  (B-C)  =D-h’  +((n/m)(B-Q) 
solving  for  h’  when  n/m  =1  (reflects  spacing  of  RWD  sensors): 

h'  =  B-2C+D 

Consequently,  the  deflection  at  that  point  on  the  pavement  due  to  the  applied  load  is: 

d  =  h-h* 

This  algorithm  is  completely  independent  of  variations  in  the  beam  angle  and  beam 
height  above  the  pavement  (remember,  h  and  h'  are  virtual  heights). 
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D  C  B  A 


D=4  C=1  B=3  A=3 


Figure  E-5  Variable  Beam  Angle  and  Height  (after  Johnson,  1995) 

Figure  E-5  provides  an  example  of  how  beam  angle  and  height  above  the  pavement  do 
not  effect  the  ability  to  measure  deflections  in  the  pavement.  In  the  top  figure  of  Figure 
E-5,  the  virtual  height,  h  can  be  calculated  as:  h  =  A-2B+C  =  5-2(2)+4  =  5 

and  the  bottom  figure  of  Figure  E-5  where  the  beam  has  moved  forward  by  one  sensor, 
the  virtual  height  h'  can  be  calculated  as:  h'  =  B-2C+D  =  3-2(l)+4  =  5 

from  this,  it  is  clear  that  the  measured  deflection,  h-h'  is  equal  to  0. 


APPENDIX  F 


TEMPERATURE  ADJUSTMENT 


F.l  INTRODUCTION 

The  temperature  adjustment  procedure  developed  as  part  of  this  thesis  is 
presented  in  an  appendix  rather  than  in  Chapter  4  with  the  load  adjustment  because  this 
procedure  is  not  required  when  using  the  regression  equations  presented  in  Chapters  3 
and  5.  The  procedure  is  presented  in  this  section  to  provide  a  record  for  future  use  in 
redeveloping  regression  equations  with  a  constant  AC  moduli  as  recommended  in 
Chapter  6. 

A  temperature  adjustment  is  required  if  the  testing  will  occur  at  other  than  a 
standard  pavement  temperature  of  25  degrees  Celsius.  Deflections  and  the  Area 
Parameter  calculated  by  those  deflections  will  increase  as  the  pavement  temperature 
increases.  The  pavement  is  stiffer  at  colder  temperatures  and  therefore  will  deflect  less 
than  at  higher  temperatures.  Section  2.6.3  in  Chapter  2  discussed  other  research  on 
temperature  effects. 

F.2  METHODOLOGY 
F.2.1  Data  Generation 

The  following  relationship  between  pavement  stiffness,  E^c,  and  temperature,  T, 
is  well  accepted  and  used  (WSDOT  Pavement  Guide,  1995): 
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Based  on  this  equation.  Table  F-1  shows  the  range  of  temperatures  associated  with  the 
asphalt  concrete  stiffness’  used  for  the  initial  120  pavement  scenarios. 


Table  F-1  Ranges  of  Eac  Used  vs  Temperature 


Eac 

Temperature 

(MPa) 

CO 

690 

58.0 

1380 

42.5 

2760 

27.0 

5520 

11.4 

11040 

O 

1 

Based  on  this  information  alone,  it  would  seem  that  temperature  differences  were  already 
included  in  the  data  that  generated  the  relationships  using  the  pavement  responses  from 
the  initial  pavement  scenarios.  However,  another  factor  affects  pavement  stiffness.  The 
level  of  fatigue  will  affect  the  pavement  stiffness.  The  temperature  vs  Eac  equation  is 
based  on  unfatigued  pavement,  however,  once  a  pavement  begins  to  fatigue  (crack)  the 
stiffness  will  also  reduce. 

In  order  to  proceed  with  generating  temperature  shift  factors,  an  assumption  was 
made  to  clarify  the  source  of  the  data.  The  original  data  from  the  120  pavement 
scenarios  is  assumed  to  be  at  a  constant  temperature  of  25  degrees  Celsius  meaning  the 
AC  stiffness  ranges  are  due  to  fatigue  levels  only.  This  assumption  clarifies  that 
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temperature  adjustments  are  not  already  embedded  in  the  data  due  to  the  range  of 
stiffness’  used. 

F.2.2  Re2ression  Analysis 

For  the  next  step,  all  data  is  assumed  to  be  for  unfatigued  pavement  at  different 
temperatures  that  correlate  to  the  AC  stiffness’  as  noted  in  Table  F-1.  In  this  manner,  the 
same  data  can  be  used  to  see  the  effects  of  temperature  on  Dq  and  Area  Parameter.  Each 
combination  of  AC  Thickness  and  Ego  (24  combinations  possible)  was  used  to  plot 
Temperature  vs  1/A  and  Temperature  vs.  Do.  These  are  the  relationships  that  correlate 
closest  to  a  straight  line  based  on  plots  on  Minitab.  The  regression  equation  for  the 
“best-fit”  line  was  generated  for  each  of  the  24  combinations  and  can  be  found  in 
Appendix  D.  The  slope  of  this  line  is  the  temperature  adjustment  for  that  combination 
and  is  shown  in  Table  F-2  for  Area  Parameter  and  in  Table  F-3  for  Dq. 


Table  F-2  Temperature  Shift  Factors  for  1/A 


Esc 

— 

— 

AC 

Thickness 

(cm) 

(MPa) 

"2 

.5 

3. 

75 

5.0 

10.0 

20.0 

40.0 

35 

X 

mifl 

fSM 

^ESM 

70 

2 

X 

wm 

[EH 

DBEH 

HQEH 

UBEM 

3 

X 

mm 

Ba 

jm 

■iiBiitia 

HQEII 

HQEH 

1  280 

3 

X 

mm 

9x 

ilM 

HQEII 

HBJbH 

HQEH 
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Table  F-3  Temperature  Shift  Factors  for  Dq 


Esg 

(MPa) 

AC  Thickness  (cm)  || 

2.5 

3.75 

5.0 

10.0 

20.0 

40,0 

35 

3.5 

5.08 

6.33 

8.29 

8.14 

6.83 

70 

2.18 

3.71 

4.83 

6.36 

6.14 

5.30 

140 

1.51 

2.93 

3.93 

5.12 

4.88 

4.35 

280 

1.22 

2.53 

3.42 

4.37 

4.13 

3.79 

F.3  RESULTS 


The  shift  factors  in  Table  F-2  and  Table  F-3  were  then  correlated  against 
subgrade  modulus.  Ego,  and  AC  thickness,  T^c,  to  see  if  a  relationship  exists  to  be  able  to 
predict  the  temperature  shift  factor  for  any  pavement  scenario.  The  best  equation  for 
each  shift  factor  follows.  Additional  equations  attempted  can  be  found  in  Appendix  D. 


logAgp  =  0.849  +  0.0361(EsJ°'-1.93 


f-1 


r2  =  92_9o/^ 


r2  =  90.8% 
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calculated,  the  adjusted  Dq  and  A  are  calculated  by  the  equations; 


Aadjusted  =  Atemp  +  (25-temp)(AsF) 
=  Do.„p  +  (25-temp)(Do,,) 
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where  Do(temp)  and  Atemp  are  the  values  measured  in  the  field  and  temp  is  the  pavement 
temperature  during  field  testing. 

The  adjusted  value  for  A  will  be  greater  that  the  measured  value  of  A  when  the 
temperature  during  testing  is  less  than  25  °  C.  The  same  relationship  will  exist  for  Dq. 


